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FOREHOIO  ■ 

Itater  vapor  strongly  Influences  the  atesphere's  ability  1»  refri^t 
electromagnetic  (EN)  waves.  In  particular*  evaporation  from  the 
surface  creates  a  vertical  gradient  of  water  vapor  that  caoi^  / 
horizontally-propagating  EN  wave  trains  to  be  refracted  downmrd.  Hhen 
an  EN  wave  train  Is  refracted  downward  with  curvature  greater  than  that 
of  the  earth's  surface,  the  wave  train  Is  said  to  be  ducted.  And  when 
ducting  of  this  kind  Is  controlled  by  evaporation,  the  zone  In  which 
ducting  occurs  Is  called  the  evaporation  duct. 

To  understand  how  EN  wave  trains  are  ducted  In  atmospheric  layers 
overlying  the  sea  surface.  It  Is  helpful  to  know  the  evaporation  duct 
height  (Z*).  the  vertical  distance  from  the  sea  surface  to  the  point  In 
the  atmosphere  where  ducting  first  ceases  (because  of  a  diminished  water 
vapor  gradient). 

To  learn  more  about  radar  ducting  (and  about  EN  wave  ducting  In 

general),  the  Ocean  Neasurements  Program  (01^)  of  the  U.  S.  Naval  - 

Oceanographic  Office  (NAVOCEANO)  contracted  with  the  Environmental 
Physics  Group  of  the  Naval  Postgraduate  School  (NPS)  to  detemlne: 

1)  some  typical  Nid-Atlantic  values.  2)  the  validity  of  a  simplified 
"bulk  method*  designed  for  routinely  monitoring  Z*  values,  and  3)  the 
parameters  most  likely  to  preswt  measurement  problems  that  could 
obstruct  the  routine  monitoring  of  evaporation  duct  characteristics. 
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I.  Sumnary 

A.  Introduction 

This  Is  a  report  on  atnospheric  measurements  made  by  the  Naval 
Postgraduate  School  Environmental  Physics  Group  aboard  the  USNS  Kane 
during  operations  In  the  North  Atlantic  In  March  of  1978.  The  primary 
goal  of  these  measurements  was  to  obtain  measurements  of  the  height  of  the 
radar  evaporative  (surface)  duct  from  surface  layer  water  vapor  and 
temperature  profiles  and  to  evaluate  the  practicality  of  using  a  simpli¬ 
fied  (the  so  called  “bulk"  method)  method  for  more  routine  estimations. 

The  bulk  method  requires  measurements  of  only  four  quantities:  the  sea 
surface  temperature  and  the  air  tenperature,  wind  velocity  and  relative 
humidity  at  some  reference  height  (usually  10  meters).  It  also  has  the 
advantage  of  considerably  less  stringent  measurement  accuracy  require¬ 
ments.  Due  to  the  extreme  difficulty  In  making  accurate  profile  mea¬ 
surements  at  sea,  the  bulk  method  Is  essentially  the  standard  for  evap¬ 
orative  duct  height  calculations. 

B.  Description  of  the  Problem 

In  a  dielectric  medium,  electromagnetic  wave  propagation  1s 
Influenced  by  the  refractive  Index.  Just  as  an  abrupt  change  In  refrac¬ 
tive  Index  causes  EM  waves  to  be  "bent",  so  does  a  continuous  gradient  In 
refractive  Index.  Since  water  vapor  Is  a  major  contributor  to  N,  there 
usually  Is  a  strong  gradient  In  N  over  the  sea  surface  because  of  evap¬ 
oration.  This  gradient  causes  EM  waves  to  be  refracted  downward.  If  a 
horizontally  propagating  wave  experiences  a  refraction  greater  than  the 
earth's  curvature,  the  wave  Is  said  to  be  ducted.  Since  the  surface  based 
duct  Is  caused  by  evaporation.  It  Is  often  called  evaporation  duct.  The 
water  vapor  gradient  (and  therefore  the  N  gradient)  over  the  ocean  typi¬ 
cally  decreases  with  height  In  the  surface  layer.  As  a  result,  at  some 
height,  the  N  gradient  has  decreased  below  the  value  required  for  ducting, 
this  Is  called  the  evaporation  duct  height,  Z*. 

Since  ducting  can  lead  to  signal  enhancements  of  many  orders  of 
magnitude,  knowledge  of  N  and  In  the  surface  Is  very  useful.  Fortunately, 

the  physics  of  the  surface  layer  Is  well  understood  and  the  specifications 
of  the  N  profile  and  Z«  can  be  reduced  to  the  determination  of  only  four 
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scaling  paramters.  This  theoretical  framework,  known  as  Monin-Obukhov 
Similarity  (MOS),  Is  based  upon  the  relation  of  the  surface  layer  profiles 
of  temperature,  wind  speed,  water  vapor  and  turbulence  with  the  surface 
fluxes  of  momentum  and  sensible  and  latent  heat.  As  one  might  expect,  the 
problem  Is  simple  In  principle  but  difficult  In  practice,  primarily  because 
the  fluxes  cannot  be  easily  measured  from  an  unstable  platform.  At  present, 
the  best  technique  Is  to  measure  the  profiles  at  several  different  heights, 
solve  the  MOS  expressions  “backwards"  to  determine  the  scaling  parameters, 
thus  determining  the  N  profile  and  Z«.  This  process,  equivalent  to  Infer¬ 
ring  the  surface  fluxes  from  the  profiles.  Is  referred  to  as  profile 
method. 

Multi-level  measurements  necessary  for  the  profile  method  are 
not  trivial.  For  example,  from  measurements  at  heights  of  10  and  25 
meters  above  the  sea  surface,  one  can  typically  expect  differences  on  the 
order  of  10%  In  wind  speed,  3%  in  relative  humidity,  and  0.2C  In  tempera¬ 
ture.  To  determine  the  scaling  parameters  to  a  10%  accuracy,  measurement 
accuracies  required  are:  wind  speed:  1%,  relative  humidity:  0.3%  and 
temperature:  0.02C.  Note  that  the  distortions  caused  by  the  presence  of 
the  ship  could  exceed  these  accuracies.  Atmospheric  measurements  even 
approaching  this  precision  require  highly  trained  scientific  personnel 
with  state-of-the-art  Instrumentation.  These  extremely  taxing  measurement 
accuracies  have  led  to  the  Interest  In  the  bulk  method. 

If  one  considers  the  MOS  expressions  for  the  gradient  as  mathe¬ 
matical  derivatives,  then  the  Integral  of  the  gradient  from  the  surface 
CZ»Zg*0)  to  some  reference  height  (Z)  specifies  the  variable  (say  tem¬ 
perature,  T(Z)).  Note  that  we  have  Introduced  two  Integration  boundary 
conditions  T(0)  and  Z^,  where  T(0)  is  the  sea  surface  temperature  and  Z^ 

Is  the  roughness  length.  If  we  assume  that  the  wind  speed  at  the  surface 
Is  zero  (U(0)  »  0)  and  the  relative  humidity  at  the  surface  is  100%  (H(0) 

»  100),  only  four  quantities  need  to  be  measured:  T(0),  T(Z),  U(Z),  and 
H(Z).  However,  an  unknown  quantity  Z^,  which  we  have  also  Introduced, 
must  be  empirically  determined. 

Average  temperatures  measured  during  the  cruise  were 
T(0)  ■  16. 2C  and  T(Z)  ■  14. 7C,  a  difference  of  1.5C  -  an  order  of  mag¬ 
nitude  greater  than  the  two  level  example  given  previously.  In  short,  the 


bulk  method  has  two  advantages:  fewer  measurements  are  required  and  an 
order  of  magnitude  greater  toleration  to  error.  This  simplicity  was 
obtained  by  Introducing  the  quantity  which  has  been  empirically 
studied  by  statistically  averaging  vast  numbers  of  over-water  measurements 
of  bulk,  profile,  and  turbulence  data.  One  of  the  primary  questions  about 
the  bulk  method  Is  the  applicability  of  an  "average"  to  a  specific 
situation.  An  even  more  serious  question  deals  with  the  validity  of 
extrapolating  the  MOS  gradient  expressions  (essentially  Z"^)  down  to  the 
surface. 

C.  Summary  of  Measurement  Results 

The  total  data  base  for  the  cruise  consists  of  360  periods  (half 

hour  averages)  of  which  40  periods  contained  near  surface  (Z»2.8m)  bouy 

profile  measurements.  A  tabulation  of  the  EM  propagation  data  (Z«,  Fm, 
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N,  dN/dZ  and  )  Is  given  In  Appendix  C.  There  are  several  ways  of 
looking  at  the  comparison  of  the  profile  or  turbulence  values  with  the 
bulk  values.  For  Instance,  the  bouy  profile  values  for  the  water  vapor 
scaling  parameter,  q«,  agreed  very  well  with  the  bulk  values,  on  the 
average,  but  the  variance  of  Individual  values  was  about  50%.  Since  a 
considerable  portion  of  this  50%  Is  due  to  uncertainty  In  the  profile 
measurements,  these  results  are  probably  rather  conservative.  The  com¬ 
parison  Is  summarized  In  the  following  table: 


QUANTITY 

SCALING  PARAMETER 

AVERAGE  AGREEMENT 

RMS  VARIANCE 

%VARIANCE 

Water  Vapor 

q* 

excellent 

.12  gm/kg 

50% 

Temperature 

T* 

poor 

.045C 

50% 

Wind  Speed 

u* 

fair 

.06  m/sec 

15% 

r  'ct  Height 

z* 

fair 

10  M 

100% 

0.  Conclusions 

Duct  heights  calculated  on  the  basis  of  bulk  measurements  are 
probably  within  a  factor  of  2  of  the  correct  value,  and  possibly  much 
closer.  The  actual  calculation  of  Z*  from  T(0),  U(Z)  and  H(Z)  can  be 
done  on  a  hand  held  programmable  calculator.  Given  the  sensitivity 
analysis  of  Section  V,  we  feel  that  the  real  problem  in  routine  appli¬ 
cations  Is  the  temperature  measurements.  The  wind  speed  can  be  obtained 
with  sufficient  accuracy  using  "eyeball"  estimates  and  a  Beaufort  scale. 


The  relative  humidity  measurements  can  probably  be  made  with  a  sling 
psychometer  from  a  well  exposed  portion  of  the  ship.  However,  there  is 
no  escaping  the  fact  that  the  air  sea  surface  temperature  difference 
must  be  measured  to  better  than  0.5C.  This  requires  well  calibrated  and 
exposed  temperature  sensors  (with  low  drift  characteristics).  The  air 
temperature  measurement  needs  to  be  made  with  a  good  aspirator  in  a 
location  reasonably  free  of  ship  influence. 

Typical  bulk  values  estimated  during  the  cruise  were  around 
15  meters.  Under  high  relative  humidity  conditions  (rain  and  fog)  Z* 
was  as  low  as  1  meter.  The  maximum  Z*  was  around  25  meters.  Consid¬ 
erable  time  variations  were  observed,  mostly  associated  with  frontal 
passages.  In  view  of  this,  measurements  should  probably  be  made  several 
times  a  day. 

One  result  of  the  theoretical  analysis  was  the  astounding 
volatility  of  Z*  under  stable  conditions  (air  warmer  than  water).  For 
the  example  with  T(Z)-T(0)=1 .4C,  a  relative  humidity  of  90%  gave  Z*  =  5 
meters  whereas  a  relative  humidity  of  75%  gave  Z*  =  60  meters.  Unstable 
conditions  showed  much  less  variability.  In  fact,  since  the  relative 
humidity  over  the  ocean  is  rarely  less  than  40%,  we  expect  Z*  to  almost 
always  be  less  than  30  meters  during  unstable  conditions.  Therefore, 
stable  conditions  produce  unusually  strong  surface  ducting  and  the  most 
variable  ducting  conditions.  Although  expected,  the  effect  is  larger 
than  we  would  have  anticipated. 

E.  General  Comments  and  Recommendations 

1.  We  feel  another  cruise  is  called  for  to  clear  up  the 
uncertainties  of  the  Kane  data.  We  are  now  "on  line"  with  improved 
temperature  and  humidity  profile  instruments,  an  IR  sea  surface  tem¬ 
perature  sensor  and  an  improved  computerized  data  logger  and  analyzer. 
One  more  cruise  could  do  a  great  deal  in  establishing  the  bulk  method 
procedures. 

2.  An  EM  propagation  analysis  of  other  NPS  data  could  be 
compared  to  the  Kane  data  or  to  climatological  predictions.  We  have  on 
file  data  from  the  Pacific,  the  North  Atlantic  and  the  Mediterranean. 

3.  More  work  needs  to  be  done  on  the  drag  coefficients 
(related  to  Z^).  The  wind  speed  coefficient  is  fairly  well  known  but 


real  knowledge  about  the  temperature  and  water  vapor  coefficients  Is 
quite  rare. 

4.  Development  of  a  simple  system  for  routine  shipboard 
measurements  of  sufficient  quality  to  allow  bulk  calculations  should 
probably  be  undertaken. 

5.  A  great  deal  of  boundary  layer  work  needs  to  be  done 
(above  the  surface  layer)  to  allow  extension  of  similarity  scaling  to 
greater  heights.  If  successful.  It  may  be  possible  to  make  estimates  of 
elevated  ducts  from  surface  and  remote  sensing  data. 
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II.  Instrumentation 

A.  Sensors  and  Placement 

Rather  than  embark  upon  a  detailed  description  of  the  NPS 
Instrumentation,  we  will  include  as  Appendix  A  a  report  (Houlihan  et  al , 
1977)  on  the  subject  and  give  only  a  brief  sketch.  The  mean  measure¬ 
ments  are  summarized  in  the  following  table: 

MEASUREMENT  INSTRUMENT  ACCURACY 


Temperature,  T 

Quartz  oscillator 

.OIC 

Relative  Humidity,  RH 

LiCl  cell 

3% 

Wind  Velocity,  V»U 

Cup  anemometer 

n 

The  relative  humidity  and  temperature  sensors  are  deployed  in 

aspirated  shelters  to  reduce  errors  from  solar  and  IR  radiation.  These 

aspirators  represent  a  weak  point  in  the  temperature  profile  measurements 

and  do  not  produce  reliable  temperature  profile  data  during  the  day. 

The  turbulence  data  was  based  upon  measurement  of  fluctuations  of  wind 

velocity  (yielding  the  rate  of  dissipation  of  turbulent  kinetic  energy, 

e)  and  temperature  fluctuations  (yielding  the  temperature  structure 

function  parameter,  Cj  ).  Constant  temperature  hot  film  sensors  were 

used  to  measure  e  and  microthermal  platinum  wires  were  used  to  measure 
2 

Cj  .  The  frequency  response  of  these  measurements  is  greater  than 
500HZ. 

The  instruments  are  deployed  at  discrete  levels:  two  levels 
on  the  NPS  tower,  one  on  the  bow  davit  and  (circumstances  permitting) 
one  on  the  buoy  (see  Figure  la).  The  following  are  the  level  heights 
(Z)  above  the  water: 

Level  #  1234 

Height,  meters  2.8  9.1  18.6  24.7 


Z2»  9.1m 
Z3-  laem 
Z4  ■24.7m 


Figure  la.  Mounting  Arrangement  on  USNS  KANE  for  NFS  Equipment. 

B.  BUOY 

The  buoy,  which  is  lowered  into  the  water  by  a  crane  on  the 
bow  of  the  ship,  is  used  to  obtain  a  measurement  very  close  to  the 
surface  (Figure  lb).  Since  the  gradients  we  wish  to  measure  vary 
roughly  as  Z"\  they  are  easier  to  measure  near  the  surface.  Unfortu 
nately,  the  buoy  can  only  be  used  when  a)  the  ship  is  stationary,  and 
b)  the  seas  are  moderate. 


III.  Theory 

A.  EM  Propagation 

1.  Refractivity  N  and  Ducting 

The  refractivity  N  (in  N  units)  of  the  atmosphere  for  EM 
wave  propagation  is 


N  =  77.6  P/T  +  3.73  x  10^  e/T^ 


(1) 


Where  P  is  the  atmospheric  pressure,  T  the  absolute  temperature  and  e  is 
the  partial  pressure  of  water  vapor.  The  vertical  gradient  of  N  is 


dN 
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^1  3Z  ^  ^2  3Z  ^  ^3 


(2) 


where  q  is  the  water  vapor  mixing  ratio  (q  *  .625  e)  in  gm/kg  and 


'1 


.3,  C2  -  7.2  and  -  -1.3.  A  negative  gradient  in  refractivity  at 


height  7.  causes  a  horizontally  propagating  EM  wave  to  be  bent  towards  the 
surface.  If  the  gradient  is  large  enough,  the  amount  of  bending  will  ex¬ 
ceed  the  curvature  of  the  earth  and  the  wave  will  be  ducted  at  the  height 
Z.  Since  the  magnitude  of  the  gradient  decreases  with  increasing  height 
in  the  surface  layer,  at  some  height  (Z*)  the  gradient  will  no  longer  be 
strong  enough  for  ducting.  Waves  propagating  below  Z«  are  ducted,  those 
above  are  not.  This  height  is  called  the  evaporative  duct  height,  and 
is  defined  by  the  critical  gradient  necessary  for  EM  wave  trapping. 


.  1C7  /,> 

37  -  -*^57  m  (3) 

Using  Eq.  21  and  noting  that  dP/dZ  »  -pg  (the  normal  hydrostatic  balance 
for  the  atmosphere),  the  is  defined  by  the  following  equation 

-.157  »  -.032  +  Cg  U  +  C3  U  (4) 

Since  3q/aZ  and  aT/aZ  are  height  dependent  (approximately  'v.Z”^  there 
will  exist  a  value  of  Z  »  Z«  such  that  Eq.  4  is  satisfied. 


T 


2.  Structure  Function  Cj^ 

Turbulence  In  the  atmosphere  produces  fluctuations  in  N  on 
a  broad  range  of  size  scales  tmi 1 imeters  to  kilometers).  These  fluctua¬ 
tions  scatter  the  EM  waves  producing  a  loss  of  resolution  and  coherence. 
If  N(x)  and  N(x+d)  represent  values  of  N  at  two  points,  separated  a  dis¬ 
tance  d,  then 

-  <(N(x)  -  >l(x+d))^>d"2/3  (5) 

For  EM  Waves 


T  o 


where  A  and  B  are  constants  is  the  water  vapor  structure  function 

and  Cjq  is  the  covariance  function.  Neglecting  the  Cj^  term  (at  sea  level) 

C  2  «  [  (77.6  X  10~^)(7.67  x  10^)j2 
N 

B.  Mon in-Obukhov  Similarity  Surface  Layer  Equations 

The  vertical  gradients  of  the  wind  speed  (U),  potential  tem¬ 
perature  (T)  and  water  vapor  mixing  ratio  (q)  are  given  by 

■  4^  *x  (5>  <8> 

Where  X  *  U,  T  or  q,  k  =  Von  Karmon's  constant  (35),  is  the  diffusi- 
vity  constant  (a^j*1.0),  X*  is  the  scaling  parameter  and  1^(5)  is  the 
MOS  stability  correction  function  (Appendix  B).  Also 

C  »  Z/L  (9) 

where  L  is  the  Monin-Obukhov  stability  length  defined  as 

“  Kg(T*^*.l8  qj 

If  we  integrate  Eq.  8  from  Z  =  to  Z,  we  obtain 

Where  i|/„(5)  is  the  profile  stability  function  (Appendix  B). 

X  2 

The  structure  function  parameter  is  related  to  X*  by  the 
function  f„(c)  (Appendix  B) 

,  xj  Z"^''^  (C)  (12) 

A  similar  form  exists  for  e„,  the  rate  of  dissipation  of  fluctuations 
of  X 

('3) 

Where  is  the  Corrsin  constant  and  E  (c)  is  the  MOS  dissipation  func- 

A  A 

tion.  In  this  format,  »  e. 
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C.  MOS  Parameter  Determination  Methods 
1.  Profile  Method 

The  profile  method  is  based  upon  Eq.  11,  expressed  in  the 

following  form 


X(Z)  =  ^  Z  -  ♦x(5)]  ♦  [X(0)  -  ^  In  Z„,] 


(14) 


This  can  be  expressed  as  a  linear  regression  of  the  form 

X(Z)  =  mY  +  b  (15) 

Where  Y  =  In  Z  'l»j{(c).  If  we  have  two  or  more  levels  of  X(Z),  we  can  fit 
the  data  to  Eq.  16,  the  slope  m,  of  the  fit  gives  X*  =  Oj^icm.  Note  that 
it  is  not  necessary  to  know  the  constant  factor,  b,  (related  to  Z^^^) 
to  find  X*.  The  actual  process  is  done  iteratively  in  the  following 
manner: 


a.  Assume  a  value  for  L 

b.  For  each  data  point  at  height  Z,  calculate  c  and 

♦xU). 

c.  Least  squares  fit  the  multi-level  data  to  obtain 

U*,  T,,  and  q^. 

> 

d.  Calculate  a  new  value  of  L  from  Eq.  10 

e.  Return  to  b)  until  L  converges 

Since  the  wind  speed  profiles  from  the  Kane  were  consi¬ 
derably  influenced  by  the  ship,  we  used  the  bulk  method  to  calculate 

u*. 

2.  Bulk  Method 


The  bulk  method  is  also  based  on  Eq.  15.  If  one  knows 
values  for  then  X*  can  be  calculated  simply  from  X(Z)  and  X(0). 
Rewriting  Eq.  15 


CX(Z)  -  X(0)) 


x’'  (In  Z/Z 


ox 


-  ^r?y) 


(16) 


We  see  the  form  of  the  standard  drag  coefficient  equation 

X*  «  (X(Z)  -  X(0)) 

Where  Cj^  Is  the  drag  coefficient  for  X 

Cv 


“XN 


[1  -  (<..<)■'  c 


XN 


ITT 


.,(01- 


(17) 

(18) 
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In  this  form,  is  related  to  the  neutral  drag  coefficient  by 


'XN 


1/2 


Ct  1C 

X 


(19) 


The  MOS  stability  parameter,  Ct  can  be  obtained  from 
Eq.  10  and  Eq.  18. 

5  .  (1  - 

°  (1  -  (o^)-’  *^(0) 


Where 


(20) 


1/2 


5  .£92.  _ C(T(Z)  -  T(0»  *  .18(q(Z)  -  a(0))1  ,3,, 

9  '  ’•UN  U(Z)^ 


In  this  case,  the  process  is  somewhat  simpler,  one  cal¬ 
culates  Cq  from  the  data  and  solves  Eq.  21  iteratively  for  c.  Given  5. 
the  MOS  parameters  are  calculated  straightforwardly  from  Eqs.  18  and  19. 
3.  Turbulence  Method 

Given  a  reasonable  estimate  of  5  (from  either  of  the  two 

previous  methods)  one  can  obtain  an  independent  value  of  the  scaling 

2 

parameter  from  turbulence  measurements  of  Cy 
(Eq.  12) 


Cq  or  e.  For  instance 


z^/^/f^(c) 


(22) 

(23) 


or  (Eq.  13) 

0/  -  eZ/Ey(c) 

D.  EM  Propagation  Equations  in  MOS  Form 

Near  sea  level,  the  gradient  of  N  can  be  written  (Eqs,  2  and  3) 

♦j(C) 


dN/dZ  »  -.032  + 


XjkZ 


(7.2  q*  -  1.3  T*) 


(24) 


The  duct  height,  Z*,  is  obtained  from  Eq.  24  by  setting 
dN/dZ  »  -.157,  the  critical  gradient  for  ducting. 

-(7.2  q*  -  1.3  T*) 
oyK  (.i2'5) 


Z*  » 


(25) 


Given  q*,  T*  and  L,  we  solve  Eq.  25  teratively  to  obtain  Z*.  Note 
that  no  ducting  will  occur  unless 

7.2  q*  -  1.3  T*  <0  (26) 
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The  refractive  index  structure  paraneter  is  calculated  from  Eqs.  7 
and  12 

C^2.^177.6»10-V67JS..V0^)]^  q.2  (27) 

IV.  RESULTS 

A.  Discussion 

Before  presenting  the  results,  it  is  worthwhile  to  consider 
several  points  about  the  absolute  validity  of  the  entire  process.  The 
various  stability  correction  functions  -  ♦(€),  f(5),  E(c)  -  are  empiri¬ 
cally  determined  functions.  In  fact,  they  are  based  on  overland  mea¬ 
surements.  The  NPS  group  has  spent  a  great  deal  of  effort  in  evaluating 
the  validity  of  these  functions  over  water  and,  so  far,  they  appear 
to  be  quite  good.  There  are  other  problems  with  some  of  the  empirical 
constants  and  the  drag  coefficients.  The  profile  method  is  dependent 
upon  the  product  oi.j.tc,  both  measured  empirically.  The  estimates  we 
have  used  are  =  1.35  and  tc  =  .35,  from  Businger  et  al  (1971),  The 
actual  value  of  both  constants  is  still  a  subject  of  controversy. 

The  bulk  method  is  dependent  upon  the  drag  coefficients 
Cjyj  and  0^^^.  For  our  analysis,  we  assumed  that  water  vapor  has 
the  same  characteristics  as  temperature,  that  is 

“T  *  ®q  •  hi  ’  ^Oq  ’  ♦t^^)  ’ 

There  are  theoretical  reasons  to  believe  that  this  is  a  good  approxi¬ 
mation  but  not  exactly  correct.  Under  these  assumptions,  Cj^  *  C^i^. 

For  the  actual  form  of  the  drag  coefficients,  we  assumed  had  the 
wind  speed  dependence  given  by  Kondo  (1975).  Based  upon  measure- 

ments  summarized  by  Liu  (1977)  and  NPS  measurements  of  Ct  ,  we  assigned 

-3 

a  value  of  »  l.l  x  10  and  solved  Tor  5  as  a  function  of  5q. 

Thus,  for  each  data  value  of  5^  we  directly  calculated  5  from  the  for¬ 
mula 

C  -  5,(1  ♦  .145,-*)  (28a) 

5  -  5„(1  ♦  .135,-*  *  .265o^-')  5,>0  (281,) 
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B.  EM  Propagation  Results 

A  listing  of  the  bulk  calculations  of  the  EM  data  is  given  in 

Table  C^l  (Appendix  Cl.  The  data  includes  calculations  of  the  Z=10 

2 

meter  values  of  N,  dN/dZ  and  C|^  as  well  as  Z*  and  Fm,  Fm  is  the  mini¬ 
mum  frequency  trapped  by  the  duct  as  is  calculated  from 

Fm  =  3.6  X  10^  (29) 

Where  Fm  is  in  GHz. 

The  following  are  two  examples  to  illustrate  the  general  be¬ 
havior  of  Z*.  Figure  2  shows  Z*  for  a  three  day  period  at  anchor  ap¬ 
proximately  80  miles  west  of  the  Straits  of  Gibraltar.  The  steady  de¬ 
crease  of  Z*  for  this  first  part  of  the  record  was  associated  with  high 
relative  humidity  and  rain.  The  rapid  increase  of  Z*  around  0200  hours 
on  3/3  corresponded  to  the  passage  of  the  front  as  indicated  by  clear¬ 
ing  skies  and  a  shift  of  the  wind  from  210°  to  260°.  Figure  3  shows 
a  very  similar  event  on  station  in  the  survey  area.  In  this  case,  the 
weather  system  was  not  strong  enough  to  produce  rain  but  fog  was  ob¬ 
served  at  1800  hours  on  3/10.  The  wind  shifted  from  200°  to  300*^  be¬ 
tween  1900  hours  on  3/10  and  0700  hours  on  3/11. 


GMT(hrs) 


Figure  2.  Z*  values  for  period  3/1  to  3/4  when  KANE  was  anchored  80  miles 
west  of  Straits  of  Gibraltar. 


24.0 


GMT(hrs) 


Figure  3.  Z*  values  for  period  3/10  to  3/12  when  KANE  was  in  survey  area. 

C.  Evaluation  of  the  Bulk  Method 
1.  Scaling  Parameters 

In  order  to  make  this  evaluation  more  meaningful,  we  have 
analyzed  the  scaling  parameters  under  those  conditions  of  maximum  profile 
accuracy.  For  instance,  the  humidity  measurement  were  not  accurate  enough 
to  give  reliable  profiles  except  when  the  buoy  was  deployed.  The  tempera¬ 
ture  profiles  were  only  analyzed  for  the  night  time  periods.  The  turbu¬ 
lence  measurements  of  U*  were  valid  during  all  periods.  A  complete  tabu¬ 
lation  of  comparisons  for  all  periods  is  given  in  Table  C-2  (Appendix  C) 
but  keep  in  mind  that  the  graphs  shown  below  are  based  on  data  selected 
from  the  table.  Incidently,  the  buoy  data  comparison  has  been  extracted 
for  condensed  viewing  and  is  shown  in  Table  C-3  (Appendix  C). 

The  comparison  of  bulk  and  profile  calculations  of  q*  are 
shown  in  Figure  4  from  the  buoy  data.  On  average,  the  comparison  is 
very  good.  The  individual  point  scatter  is  about  50%.  For  this  graph, 
and  those  following,  the  error  bars  are  for  the  mean  estimates  and  tne 
number  represents  the  number  of  data  points  in  each  bin. 
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-a4  -<X3  -aa  -ai  -ao 

BULK  q«(gm/kg) 


Figure  4.  Calculation  of  water  vapor  mixing  ratio  scaling  parameter, 
profile  determination  vs.  bulk,  bouy  data  only. 

The  comparison  of  bulk  and  profile  calculations  of  T*  are  shown 
In  Figure  5.  These  results  are  very  poor  Indeed.  It  Is  not  clear  If 
the  lack  of  correlation  Is  due  to  measurement  problems  sr  Is  In  fact 
failure  of  the  bulk  method. 


•-  -0.1 
UJ 


S  -0.2 


-0.4 


■0.2  -0.1  0.0 
BULK  T*  CO 


Figure  5.  Calculation  of  potential  temperature  scaling  parameter,  T*, 
profile  determination  vs.  bulk,  night  time  data. 

The  comparison  of  the  bulk  and  turbulence  calculations  of  U*  are 
shown  in  Figure  6.  Although  the  correlation  is  not  perfect,  it  is  pro¬ 
bably  good  enough  for  bulk  determinations  of  Z*.  The  complicated  velo¬ 
city  dependence  of  that  was  used  for  these  U*  bulk  values  did  not 
give  significantly  better  comparison  than  assuming  a  velocity  indepen- 
dent  value  of  ®  1.3  x  10  . 
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Figure  6.  Calculation  of  friction  velocity,  U*,  turbulence  determin*- 
tion  vs.  bulk,  all  data, 

2.  Duct  Height 

Evaporation  duct  height,  Z*,  was  calculated  using  the 
bulk  and  profile  methods  for  all  data.  A  comparison  of  these  results 
is  given  in  Figure  7. 
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Figure  7.  Calculation  of  evaporation  duct  height,  Z*,  profile  deter 
mination  vs.  bulk.  The  open  circles  represent  all  data,  the  circled 
x's  are  the  buoy  data  only. 

Although  the  correlation  is  only  fair,  it  should  be  noted 
that  the  bins  containing  most  of  the  points  (187  out  of  232)  fall 
within  the  expected  mean  error.  Also  encouraging  is  the  much  better 
agreement  obtained  from  the  buoy  data. 

V.  SENSITIVITY  ANALYSIS 
A,  Data  Distribution 

Rather  than  consider  many  possible  combinations  of  the  me¬ 
teorological  parameters,  the  sensitivity  analysis  will  be  done  with¬ 
in  the  framework  of  typical  values  measured  during  the  cruise.  A 
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sumnary  of  the  overall  average  data  Is  given  below. 


QUANTITY 

SYMBOL 

AVERAGE 

a_ 

Sea  Surface  Temp. 

TS 

16.25°  C 

0.40®  C 

Air  Temperature 

T(10) 

14.67°  C 

1.19°  C 

Relative  Humidity 

H(10) 

70% 

17% 

Wind  Speed 

u(io) 

7.57  m/sec 

2.84  m/sec 

Mixing  Ratio 

q(10) 

7.29  gm/kg 

1.7  gm/kg 

B.  Sensitivity  of  Z* 

The  sensitivity  analysis  was  performed  usjng  the  meteorologi 
cal  parameters  as  inputs  into  the  bulk  expressions.  The  predictions 
are  shown  in  the  form  of  Z*  as  a  function  of  relative  humidity.  Of 
course,  one  should  keep  in  mind  the  extreme  rarity  of  relative  humidi¬ 
ties  (at  Z=10  meters)  less  than  40%  over  the  ocean.  On  each  of  the 
graphs  shown,  the  solid  line  represents  the  "average"  curve  (TS-16.2, 
T(10)-14.7  and  U(10)=’7.6  m/sec).  Since  the  average  relative  humidity 
was  70%,  we  expect  an  average  Z*  =  16  meters.  Figure  8a  shows  the  ef¬ 
fects  of  velocity  variations  on  the  order  of  the  standard  deviation 
found  during  the  cruise.  The  effect  of  sea  surface  variation  (at  con¬ 
stant  air-sea  temperature  difference)  is  given  in  Figure  8b.  Z*  is 
moderately  insentive  to  both  parameters.  The  corresponding  graph  for 
sensitivity  to  the  air-sea  temperature  difference  is  noticeably  differ 
ent  (Figure  9).  Not  only  can  very  large  values  of  Z*  occur  for  stable 
conditions  (AT>0)  but  note  the  very  steep  variability  of  Z*  between 
60%  and  90%  humidity.  Of  course,  one  must  bear  in  mind  that  these  are 
surface  layer  extrapolations,  the  maximum  height  of  validity  is  proba¬ 
bly  less  than  50  meters.  Nevertheless,  these  results  are  quite  sig¬ 
nificant. 


Figure  8.  Duct  height  vs.  relative  humidity  for  typical  meteorological 
conditions  found  during  the  Kane  cruise;  a)  sensitivity  to  wind  speed 
b)  sensitivity  to  sea  surface  temperature. 


Figure  9.  Duct  height  vs.  relative  humidity  for  typical  meteorological 
conditions  found  during  the  Kane  cruise:  sensitivity  to  the  air-sea 
temperature  difference .(AT  ■  T(Z)  -  T{0)). 


The  large  difference  in  stable  and  unstable  atmospheric 
evaporation  ducting  properties  is  a  direct  result  of  the  behavior  of 
the  MOS  gradient  equations.  Let  us  assume  that  very  near  the  surface 
we  have  an  N  gradient  great  enough  to  permit  ducting.  Since  the  gra¬ 
dient  very  near  the  surfaces  decreases  with  height  as  Z"\  we  expect 
the  gradient  to  decrease  with  increasing  height  until  Z=Z*  and  the  cri¬ 
tical  gradient  is  reached.  The  atmospheric  stability  enters  into  the 
problem  through  its  influence  on  the  height  dependence  of  the  gradients. 

Under  unstable  conditions,  the  convective  mixing  quickly  removes  the 

-3/2 

gradients  as  Z  increases,  leading  to  a  z  '  gradient  height  dependence. 
Consequently,  the  critical  gradient  is  reached  at  a  relatively  low  Z. 
Under  stable  conditions,  the  bouyancy  forces  are  retarding  the  turbu¬ 
lent  mixing  as  Z  increases,  leading  to  gradients  that  approach  a  con¬ 
stant  independent  of  height.  Therefore,  the  gradient  decreases  more 
slowly  under  stable  conditions  and  the  critical  gradient  can  be  reached 
at  much  greater  heights. 

VI,  CONCLUSIONS 

A.  Variability  of  Z* 

Typical  values  of  Z*  obtained  during  the  cruise  ranged  from 
around  2.0  meters  to  25  meters  with  a  mean  value  of  15  meters.  Trans¬ 
missions  at  10  GHz  would  have  been  ducted  approximately  70%  of  the 
time,  20  GHz  approximately  85%  of  the  time.  The  two  factors  having 
the  greatest  influence  on  Z*  are  the  relative  humidity  and  the  air-sea 
temperature  difference.  In  one  instance  the  duct  height  changed  from 
3.0  meters  to  20  meters  in  about  four  hours  following  a  frontal  passage. 

B.  Validity  of  the  Bulk  Method 

The  profile  and  turbulence  data  from  this  cruise  do  not  give 
a  clear  validation  of  the  bulk  method.  Although  the  water  vapor  and 
wind  speed  scaling  parameters  were  reasonably  well  correlated  in  the 
mean,  there  was  considerable  scatter  (on  the  order  of  50%  and  15%  res¬ 
pectively)  for  individual  measurements.  Comparisons  for  T*  and  Z*  were 
definitely  a  disappointment.  We  feel  that  much  of  the  discrepancy  can 
be  attributed  to  inaccurate  temperature  and  humidity  profiles.  Since 
both  of  these  systems  have  been  upgraded  since  the  Kane  cruise,  another 
cruise  should  be  made.  Based  upon  comparison  of  turbulence  T*  (from  Cj^) 
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with  bulk  T*  from  other  NPS  data,  we  feel  that  the  bulk  T*  are  much  better 
than  indicated  by  the  Kane  data. 

C.  Implications  of  the  Sensitivity  Analysis 

The  sensitivity  analysis  showed  that  the  wind  speed  and  sea  sur¬ 
face  temperature  contributions  to  Z,,,  though  significant,  were  not  highly 
critical.  The  relative  humidity  and  air-sea  temperature  difference  were 
definitely  critical  factors.  Given  the  standard  measurement  capabilities 
of  these  two  quantities,  it  is  obvious  that  the  air-sea  temperature  dif¬ 
ference  will  be  the  primary  roadblock  to  implementing  a  program  of  rou¬ 
tine  measurement  of  Z*  either  from  ship  data  or  satellite  data.  Another 
result  of  the  study  was  the  large  values  of  Z^  possible  under  stable  con¬ 
ditions  (stable  conditions  predominate  off  the  east  coast  of  the  U.S.  and 
Canada  in  the  summer). 
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ABSTRACT 


Programs  of  Investigation  of  optical  propagation,  marine  fog,  and 
boundary  layer  properties  In  the  marine  environment  have  led  to  the 
development  of  a  shipboard  system  for  measuring  atmospheric  stability 
and  turbulence.  This  research  requires  verv  accurate  multilevel 
measurements  of  mean  and  fluctuating  temperature,  wind  velocity,  and 
humidity.  This  paper  is  a  description  of  the  system  including  the 
following  aspects:  1)  instrumentation,  2)  calibration,  3)  analysis 
techniques,  and  4)  examination  of  the  data  and  evaluation  of  the 


system. 
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I.  Introduction 

There  is  currently  considerable  interest  in  curoulent  transport  pro¬ 
cesses  in  Che  marine  boundary  layer.  Two  problems  which  are  of  parti¬ 
cular  Importance  concern  some  turbulence  properties  and  their  influence 
on  optical  propagation,  and  Che  correlation  between  the  transport  of  heat 
and  moisture  and  Che  occurrence  of  marine  fog.  In  Che  optical  problem, 
turbulence  causes  beam  wander  and  spreading  and  scintillation  (twinkling), 
which  cause  energy  density  loss,  resulting  in  information  loss  in  the  case 
of  communication  systems.  The  occurrence  of  marine  fog  is  strongly  depen¬ 
dent  on  Che  molsCure/heaC  ratio  in  Che  atmosphere  and  Che  turbulent  trans¬ 
port  of  these  quantities  is  of  critical  Importance.  The  goal  of  present 
investigations  of  turbulence  properties  is  to  enable  predictions  of  Che 
relevant  atmospheric  properties  on  Che  basis  of  relatively  simple  obser¬ 
vations  Chat  are  currently  being  used  by  weather  centers,  k  sophisticated 
secies  of  experiments  is  necessary  at  this  stage  of  development  because 
of  the  need  to  develop  an  adequate  parameterization  of  the  problem  and 
because  of  the  complications  inherent  in  model  development  and  verification. 

The  Naval  Postgraduate  School  research  vessel  R/V  Acanla  is  the 
primary  platform  for  overwacer  measurements  and  it  is  the  system  which  is 
installed  on  this  ship  which  is  described  here.  A  four-level  system  has 
been  developed  to  measure  both  mean  and  fluctuating  quantities.  This 
system  is  augmented  by  a  vlslometer,  a  sea  surface  temperature  probe,  an 
acoustic  sounder,  and  at  times  by  balloona  kytoons  or  kites  to  make  ele¬ 
vated  measurements. 

Since  a  ship  is  a  non-stable  platform  there  is  always  some  question 
as  to  the  validity  of  data  obtained  when  turbulence  parameters  are  being 
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Investigated.  It  Is  not  possible  to  measure  turbulence  directly  at  low 
frequencies  (<lHz}  aboard  ship  because  the  ship  motion  interferes  in 
this  region.  At  higher  frequencies  valid  data  can  be  obtained,  and  the 
experimental  and  analysis  techniques  described  herein  exploit  this  fact. 
An  additional  problem  with  shipboard  measurements  is  the  large  physical 
size  of  a  ship  and  the  disturbance  of  the  local  air  flow  by  the  ship's 
presence.  Fortunately,  the  R/V  Acanla  Is  a  fairly  small  and  narrow  ship 
and  this  effect  is  not  large.  Moreover,  proper  placement  of  sensors  to 
take  data  only  when  the  wind  direction  is  favorable  can  minimize  inter¬ 
ference  effects  on  the  data.  In  this  paper  evidence  is  presented  to  show 
that  neither  ship  motion  nor  physical  presence  invalidate  the  data 
displayed . 


II.  Shipboard  Configuration 

The  shipboard  sensor  mounting  arrangement  is  shown  in  Figure  A-1.  The 
sensors  for  mean  and  fluctuating  wind,  temperature  and  humidity  were 
mounted  on  two  towers,  located  approximately  2  feet  and  16  feet  behind 
the  tip  of  the  bow.  The  heights  of  the  sensors  above  the  sea  surface  were: 


Level  1 


Level  2 


4.2  meters 
6.6  meters 


Level  3 


7.6  meters 


Level  4 


13.9  meters 


The  following  sensors  were  installed  at  the  levels  indicated  on  the  masts: 


Mean  Wind 
Mean  Temperature 
Mean  Humidity 
Wind  Fluctuation 
Temperature  Fluctuation 
Humidity  Fluctuation 


Cup  Anemometer  All  Levels 

Quartz  Thermometer  All  Levels 


Li  Cl  Cell 
Hot  Wire 
Cold  Wire 
Lyman  -  a 


All  Levels 
All  Levels 
All  Levels 
Levels  1  and  3 
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The  heights  and  locations  of  the  sensors  were  chosen  In  such  a  way 
to  minimize  shipboard  Influence  and  to  measure  over  a  wide  enough  height 
difference  so  that  small  gradients  In  the  parameters  would  be  detected. 

A  subsequent  improvement  that  allows  data  to  be  obtained  much  closer  to 
the  sea  surface  will  be  described  at  the  end  of  this  section. 

The  acoustic  sounder  was  mounted  on  the  fantall  of  the  ship  since 
this  Is  the  area  that  Is  most  noise  free.  The  vibrations  of  the  deck 
when  the  ship  Is  underway  produce  enough  noise  that  useful  signals  can 
only  be  obtained  up  to  about  400  meters. 

The  vlslometer  was  mounted  on  top  of  the  bridge  on  a  S-foot  high 
support.  This  position  placed  the  instmment  In  an  uncontaminated  air¬ 
flow  when  the  relative  wind  was  within  90°  of  the  bow. 

The  sea  surface  temperature  sensor  was  located  approximately  4  feet 
from  the  starboard  side  of  the  ship  and  floats  within  one  Inch  of  the 
surface.  The  starboard  side  was  chosen  because  the  ship  water  exhausts 
were  located  on  the  port  side. 

In  order  to  study  near  surface  effects  and  to  take  advantage  of 
larger  gradients  in  the  near  surface  region,  a  buoy  systos  (Figure  A-2) 
was  developed  for  deployment  off  the  bow  of  the  Acanla  (Corbin,  1977) . 

A  heavy  weight  was  suspended  by  cable  through  a  hollow  mast  in  the  buoy, 
thus  allowing  the  buoy  to  be  rigidly  controlled  In  the  horizontal  dir¬ 
ection  but  to  freely  slide  up  and  down  with  waves. 

III.  Instrumentation 

A.  Mean  System 

A  block  diagram  of  the  mean  system  Is  shown  in  Figure  A-3.  The 
system  consisted  of  two  primary  components,  viz.,  the  sensors  and  their 
associated  electronics,  and  the  data  logging  and  readout. 
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The  data  logger  was  a  mlcroprogrammable  Integrated  data  acquisition 
systan  (MIDAS)  which  was  developed  at  NPS  specifically  for  acquiring 
shipboard  meteorological  data.  MIDAS  utilized  an  Intel  Model  8008 
Central  Processor  to  control  the  sampling,  averaging  and  recording  of 
mean  meteorological  data.  All  software  programs  were  written  in  PL/M 
to  facilitate  the  writing  of  self-documei^tlng  programs  (Plunkett,  1976). 

The  operator  was  interfaced  with  the  system  via  teletype  for  full 
duplex  input/output  communications  and  program  control.  The  operator 
exercised  control  over  the  sample  start-time  and  the  number  of  samples 
to  be  averaged  before  outputting.  The  operator  could  also  alter  the 
preset  sample  list  by  adding  or  deleting  various  sensors  as  they  came 
on  line  or  become  Inoperable.  Once  initiated,  the  system  was  fully 
automated  to  sample  the  tailored  list  of  sensors  every  30  seconds  and 
periodically  printed  output  values  averaged  over  a  selected  interval 
from  one  to  sixty  minutes. 

Output  values  were  printed  on  the  teletype  in  columnlzed  format 
with  the  time  of  print  as  a  leader.  The  teletype  contained  a  paper 
tape  punch  that  could  be  activated  by  the  operator  to  produce  a  data 
copy  concurrent  with  the  printout.  A  magnetic  cassette  tape  recorder 
was  integrated  into  the  system  as  a  third  data  output  device.  This 
cassette  is  to  be  interfaced  to  a  Hewlett-Packard  Model  9831  Portable 
Computer  so  that  profile  and  gradient  flux  estimates  can  be  performed 
onboard  automatically,  using  the  BASIC  programming  capability  of  the 
HP  9831.  Presently  data  cards  are  punched  from  the  paper  tape  output 
and  processed  on  the  IBM  360  System  at  NFS. 

Mean  wind  speeds  were  measured  using  Thomthwaite  Cup  Anemometer 
units  featuring  lightweight  plastic  sensors  that  ensured  low  wind  speed 
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response.  The  pulse  outputs  (one  pulse  per  revolution)  of  the  four 
anemometer  units  were  fed  directly  to  MIDAS  where  they  were  averaged 
and  printed  on  the  teletype  at  the  conclusion  of  each  preselected  time 
period.  Revolution  counts  were  also  displayed  on  Hewlett-Packard  Model 
5221A  Counters  for  quick-look  and  calibration  checks.  The  counts  were 
converted  to  wind  speed  using  the  calibration  parameters  supplied  by 
Thomthwalte,  which  were  checked  periodically  in  a  wind  tunnel. 

The  temperature  transducers  utilized  were  Hewlett-Packard  Model  2833 
Oceanographic  (quartz  crystal)  Sensors,  with  a  resolution  of  0.01°C. 

The  sensors  and  accompanying  booster  oscillators  were  sequentially  multi¬ 
plexed  through  a  Hewlett-Packard  Model  2801  Quartz  Thermometer  Readout 
Unit  for  monitoring. 

The  humidity  sensors  utilized  were  Hygrodynamics,  Inc.,  Model  1818U 
Wide  Range  Sensors,  which  featured  an  Integral  thermistor  probe  with  a 
resolution  of  O.l^C.  These  units  were  sequentially  multiplexed  through 
a  Hygrodynamics  Digital  II  Readout  Unit  for  monitoring  both  relative 
humidity  and  temperature  data. 

The  temperature  probes  in  the  humidity  sensor  arrangements  provided 
a  check  capability  on  the  more  sensitive  quartz  thermometer  sensors 
colocated  with  the  humidity  sensors.  While  the  resolution  demanded  for 
gradient  measurements  (0. 01*^0  could  only  be  satisfied  by  the  quartz 
thermometer  system,  the  availability  of  back-up,  bulk  temperatures  from 
the  humidity  units  was  useful  on  several  occasions  where  conflicting 
profile  trends  were  discerned  In  quartz  temperature  readings. 

Temperature  and  humidity  sensor  outputs  at  each  of  the  four 
measurement  levels  were  multiplexed  to  their  respective  readout  units 
using  an  WPS  developed  level  selector  which  operated  under  the  timing 
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control  of  the  MIDAS  system.  The  level  selector  fulfilled  two  essential 
requirements.  First,  It  allowed  one  readout  unit  to  serve  the  entire 
profile  array,  thus  decreasing  the  capital  Investment  In  equipment. 
Secondly,  the  level  select  unit  contained  the  power  elements  for  each 
level  in  the  array.  Thus,  when  Individual  sensors  were  selectively 
connected  to  their  respective  readout  unit,  only  signal  levels  were 
switched  at  each  junction.  Maintaining  constant  powering  at  each  level 
rather  than  switching  In  power  concomitant  with  signal  level  switching 
kept  line  transients  to  a  minimum  and  Insured  proper  time  constant 
response  for  the  temperature  and  humidity  equipment. 

The  quartz  temperature  sensors  and  oscillator  circuit,  and  the 
humidity  sensors  at  each  level  were  contained  within  C.C.  BREIDERT  CO., 
Alr-X-Houster  Aspirated  Weather  Shelters.  The  mounting  arrangement  Is 
shown  in  Figure  A-4.  The  shelters  served  two  purposes,  viz,,  protecting 
the  sensors  from  the  marine  environment,  and  reducing  radiation  effects. 
Radiation  effects  are  especially  severe  in  this  application  since  some 
of  the  sensors  are  above  the  water  and  some  above  the  deck  of  the  ship, 
which  is  heated  by  the  sun  and  radiates  strongly.  Because  of  radiation 
from  the  deck  the  shelters  were  modified  by  placing  a  shield  at  the  base 
of  the  unit  to  reduce  radiation  from  below. 

B.  Fluctuation  System 

A  block  diagram  of  the  fluctuation  system  Is  shown  In  Figure  A-5. 
It  Is  immediately  apparent  that  there  are  several  output  devices,  allow¬ 
ing  several  methods  of  data  analysis.  The  analog  tape  recorder  was  used  to 
record  the  basic  signals,  which  could  then  be  analyzed  in  the  laboratory 
subsequent  to  a  field  trip.  The  analog-digital  converters  and  the 
averaging  and  storage  were  contained  in  the  MIDAS  system  described  above 


and  chelr  results  were  printed  on  the  teletype  at  the  same  time  as  the 
mean  system  data.  Averaging  could  also  be  accomplished  directly  on  the 
strip  chart  records. 

All  sensors  were  mounted  on  wind  vanes,  with  one  vane  per  level, 
to  maintain  correct  alignment  with  the  relative  wind  directions. 

The  sensors  for  wind  speed  fluctuations  were  Thermo-Systems  model 
1210W-T1.5  probes  (4.5  micron  tungsten)  operated  as  hot  wires  with  a  50% 
overheat.  The  power  and  bridge  for  the  probes  was  the  Thermo-Systems 
model  1054  constant  temperature  anemometer,  the  output  of  which  was 
processed  by  a  Thermo -Systems  model  1057  signal  conditioner  to  remove 
the  dc  level.  The  dc  level  could  cause  the  tape  recorder  to  saturate 
when  large  fluctuation  signals  occurred.  The  signal  conditioner  also 
contained  band  pass  filters,  which  were  set  for  a  range  of  0.2  Hz  to 
1  kHz,  reducing  electronic  noise.  The  50%  overheat  used  was  high  enough 
so  that  the  system  response  to  temperature  fluctuations  was  well  below 
Che  electronic  noise  level.  The  hot  wire  sensors  were  oriented  with  the 
axis  of  the  wire  in  Che  vertical  direction.  With  this  orientation  the 
sensor  was  insensitive  Co  Che  vertical  component  of  the  wind,  so  that 
only  horizontal  fluctuations  were  recorded. 

The  sensors  for  air  temperature  fluctuations  were  Thermo-Systems 
model  1210W-P.8  (2.5  micron  platinum)  operated  as  cold  wires.  An  ac 
resistance  bridge,  Sylvania  model  140  thermosonde,  which  operated  at 
3  kHz,  was  adapted  for  shipboard  use  and  acted  as  the  detector.  With 
this  unit  data  could  be  obtained  up  to  1  kHz.  It  is  necessary  to  use 
an  ac  bridge  as  Che  detector  so  Chat  high  sensitivity  could  be  obtained 
with  a  very  low  current  through  the  sensor.  If  the  temperature  of  the 
sensor  was  elevated  due  to  bridge  current  Che  system  would  sense  wind 


speed  fluctuations  as  well  as  temperature  fluctuations.  The  bridge  had 
a  long  time  constant,  automatic  balance  circuit  so  that  slow  drifts  of 
ambient  temperature  were  compensated. 

The  thermosonde  configuration  was  such  that  external  resistances 
would  be  used  in  both  arms  of  the  bridge.  This  allowed  one  to  use  a 
single  sensor  balanced  by  a  precision  decade  resistance,  or  to  use  two 
cold  wire  sensors  In  the  bridge.  Thus,  two  methods  were  available  for 
evaluating  temperature  fluctuations:  (1)  A  single  sensor,  the  signal 
being  analyzed  by  a  spectrum  analyser-T'  mode,  and  (2)  two  sensors  spaced 
at  a  fixed  separatlon*dT  mode, the  resulting  signal  being  processed  by 
an  SMS  module,  with  the  average  either  recorded  on  a  strip  chart  or  by 
the  MIDAS  system. 

Figure  A-5  shows  that  the  wind  speed  fluctuation  signals  were  also 
processed  by  an  SMS  module  after  the  signal  had  passed  through  a  band 
pass  filter.  The  band  pass  filter  operated  in  the  time  domain  and  per¬ 
formed  a  function  equivalent  to  the  spatial  filtering  accomplished  by 
the  probe  separation  utilized  in  the  AT  mode  for  temperature  fluctuations. 
Spatial  filtering  could  not  be  used  for  wind  speed  fluctuations  because 
It  was  not  possible  to  operate  two  hot  wires  by  a  single  bridge.  Using 
two  bridges  and  a  different  circuit  Introduced  too  much  error  into  the 
process. 

A  final  arrangement  for  evaluating  fluctuation  signals  led  to  the 
differentiation  of  the  bridge  output  prior  to  KMS  analysis.  Dif¬ 
ferentiation  necessarily  emphasizes  high  frequency  components  so  that 
high  frequency  noise  becomes  a  problem.  To  alleviate  this  problem  a 
1  kHz  active  filter  with  a  ^8  dB  rolloff  was  used  but  noise  problems 
were  still  severe  enough  so  chat  no  satisfactory  results  were  obtained 


with  this  method. 


Methods  for  analyzing  data  gathered  utilizing  the  various  experi¬ 
mental  techniques  will  be  described  In  Section  VI. 

Humidity  fluctuations  were  detected  with  a  Lyman-a  cell.  Difficul¬ 
ties  were  experienced  In  keeping  this  system  In  operation  In  the  rigors 
of  a  marine  environment  and  no  results  will  be  reported  here. 

C.  Auxiliary  Equipment 

An  Aerovironment  model  300  Acoustic  Radar  was  installed  on 
the  R/V  Acanla  In  order  to  monitor  the  height  of  the  temperature  Inver¬ 
sion.  The  sounder  transmitted  a  1600  Hz  pulse  of  sound  which  was 
scattered  by  small  scale  ambient  temperature  Inhomogeneltles .  Using  a 
facsimile  recording  device,  the  scattered  return  sound  waves  were  then 
detected  to  produce  the  observed  signal.  The  large  temperature  gradient 
and  strong  turbulence  at  the  base  of  an  inversion  produced  a  strong 
reflected  signal,  so  that  the  height  of  the  inversion  could  be  readily 
determined.  It  was  possible  to  measure  the  strength  of  thermal  turbu¬ 
lence  with  the  acoustic  sounder,  but  noise  and  calibration  problems  made 
such  measurements  impossible  with  our  installation. 

A  Meteorology  Research,  Inc.  model  1580A  Fog  Vlslometer  was  Installed 
on  the  bridge  of  the  R/V  Acanla.  The  Instrument  measured  the  amount  of 
light  scattered  from  particulates  and  had  a  usable  visibility  range  of 
8  m  to  13  km.  The  unit  could  be  rotated  so  that  the  relative  wind  was 
incident  directly  on  the  open  side  of  the  scattering  region.  Rotation 
was  necessary  when  the  sun  was  low  on  the  horizon  in  order  to  keep  sun¬ 
light  out  of  the  receiving  optics. 


39 


IV.  Calibration 


A.  Mean  System 

1.  Temperature 

The  mean  temperature  sensors  (HP  model  2833  Quartz)  were  cal¬ 
ibrated  using  standard  laboratory  thermometric  techniques.  Since  the 
system  was  used  to  measure  the  gradient  of  a  temperature  profile,  the 
relative  consistency  of  the  sensors  was  considered  more  important  than 
the  absolute  temperature  calibration.  All  five  sensors  were  greased  and 
screwed  into  a  3 -inch  aluminum  cylinder  fitted  with  a  center-mounted 
thermocouple.  The  aluminum  cylinder  was  0-rlng  fitted  into  a  copper 
pipe  long  enough  (24  inches)  to  permit  colling  of  the  majority  of  the 
armored  sensor  cable  in  the  temperature  controlled  environment.  The 
pipe  and  cylinder  unit  were  ismtersed  In  a  15-llter  glass  Dewar  flask  filled 
with  a  mixture  of  water  and  ice.  The  cylinder  and  pipe  combination 
allowed  the  sensors  and  cables  to  be  surrounded  by  the  bath  without 
actually  getting  wet.  The  large  mass  and  excellent  thermal  conductivity 
of  the  copper  and  aluminum  eliminated  transient  fluctuations  and  reduced 
Inter-sensor  temperature  differences  to  a  few  thousandths  of  a  centi¬ 
grade  degree. 

The  normal  calibration  procedure  was  to  cool  the  sensors  to  about 
15°C  while  monitoring  their  Indicated  temperatures  with  printout  from 
MIDAS  at  one  minute  Intervals.  Once  equilibrium  was  reached,  the  sensor 
oscillators  were  adjusted  so  that  (if  possible)  all  sensors  read  within 
.01°C  of  each  other.  The  temperature  was  then  varied  in  2-3°C  incre¬ 
ments  by  adding  ice  or  warm  water  and  noting  the  MIDAS  printout  values 
when  a  new  equilibrium  was  reached.  Initially,  the  temperatures  were 
compared  to  the  thermocouple  values  (data  set  1/7/76,  Table  A-1)  but 
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TABLE  A-1  Me2ui  Temperature  Calibrations 


T,  “C 


2/9/77 


TC 

0 

1 

2 

3 

4 

.01 

.02 

.00 

.02 

-.01 

.02 

5.25 

4.94 

4.93 

4.94 

4.92 

4.94 

7.50 

7.24 

7.24 

7.25 

7.23 

7.24 

9.45 

9.21 

9.22 

9.22 

9.21 

9.21 

11.95 

11.64 

11.66 

11.65 

11.65 

11.65 

13.35 

13.05 

13.08 

13.06 

13.06 

13.06 

14.85 

14.53 

14.56 

14.54 

14.55 

14.54 

23.25 

22.86 

22.89 

22.87 

22.88 

22.87 

8.61 

8.60 

8.62 

8.61 

8.62 

13.14 

13.14 

13.15 

13.15 

13.16 

20.54 

20.50 

20.51 

20.52 

20.53 

this  was  later  deemed  unnecessary.  One  problem  with  the  system  developed 
due  to  using  quartz  sensors  with  different  calibration  properties.  In 
data  set  1/7/76  (Table  A'-l)  sensors  1  and  3  are  from  a  different  batch 
than  0,  2  and  4.  This  led  to  Increasing  disagreement  at  the  extremes 
of  the  calibration  range.  This  problem  was  alleviated  by  ordering  more 
sensors  and  specifying  characteristics  to  match  sensors  0,  2  and  4  (data 
set  2/9/77,  Table  A-1).  With  this  Improvement,  the  accuracy  of  a  single 


sensor  (relative  to  the  average  of  all  five)  was  no  worse  than  ±  .02°C 
over  a  reasonable  operating  range.  Note  that  this  is  the  accuracy 
of  the  temperature  sensor  only  and  does  not  Include  errors  introduced  by 
the  aspirators.  The  relative  aspirator  temperature  error  depends  most 
strongly  on  the  amount  and  uniformity  of  solar  heating;  this  error  Is 
quite  small  at  night  or  under  heavy  clouds  but  may  be  as  great  as 
±  .1°C  In  bright  sunlight  with  light  winds. 

2.  Humidity 

To  determine  their  accuracy  and  reproducibility,  detailed  exam¬ 
ination  In  laboratory  conditions  has  been  made  of  the  calibration  char¬ 
acteristics  of  the  Hygrodynamics  Digital  Hygrosensor.  The  calibrations 
were  based  on  ambient  humidity  variations  (from  24Z  RH  to  50Z  RH  over  a 
one  month  period) ,  as  measured  with  a  Bendlx  Psychron  Wet-Dry  Bulb 
Psychrometer,  and  controlled  humidity  (from  40Z  RH  to  99Z  RH) ,  as  meas¬ 
ured  by  a  wet-dry  thermocouple  in  a  2-x  3-x  3-foot  fog  chamber  adapted  for 
these  experiments.  The  air  circulated  through  the  chamber  could  be  pre- 
humldlfled  in  a  column  of  wet  glass  beads  giving  nominal  humidity  of 
about  80Z.  High  humidities  (up  to  99Z  RH)  could  be  obtained  by  Intro¬ 
ducing  varying  amounts  of  water  in  the  form  of  a  fine  mist  Into  the 
chamber.  In  general,  the  humidity  In  the  chamber  could  be  stabilized 
for  periods  of  several  hours  and  was  uniform  throughout  the  chamber  to 
within  IZ  RH. 

The  hygrosensor  was  a  Dunmore-type  lithium  chloride  sensor  whose 
resistance  varies  in  proportion  to  the  relative  humidity  to  which  It 
is  exposed.  It  Is  a  slow-response  sensor  and  Is  generally  utilized  in 
mean  systm  measurements.  The  standard  procedure  of  calibrating  these 


sensors  In  a  small  chamber  with  a  saturated  salt  solution  of  known  vapor 
pressure  was  not  satisfactory  for  a  program  of  calibration  of  eight  sen¬ 
sors.  This  was  primarily  due  to  a  lack  of  confidence  in  our  ability  to 
reproduce  humidity  conditions  In  consecutive  sensor  calibrations.  Attempts 
to  calibrate  all  eight  sensors  simultaneously  In  a  large  volume  chamber 
were  complicated  by  tmnperature  drifts  that  made  It  difficult  to  assure 
equilibrium  conditions  necessary  to  the  accuracy  of  the  technique.  As 
a  result.  It  was  decided  to  use  the  humidity  controlled  chamber  previously 
described.  The  quantitative  results  of  a  series  of  calibrations  are  given 
in  Table  A-2;  In  general,  the  eight-sensor  average  agreed  with  the  psychro- 
meter  standard  to  -0.4  ±  2.9X  RH  for  measurements  taken  with  equilibrium 
times  of  order  one  hour.  The  individual  total  humidity  accuracy  of  the 
sensor  was  of  less  Importance  than  the  relative  consistency  within  the 
group  of  sensors,  since  they  were  used  in  a  four-level  system  to  determine 
humidity  gradients.  In  this  respect,  they  could  be  calibrated  in  the 
laboratory  to  t  IZ  RH,  subject  to  assurances  of  sufficient  time  for 
response. 

The  suitability  of  the  hygrosensors  for  shipboard  measurement  of 
humidity  gradients  was  reduced  by  two  Inherent  properties  of  the  device. 

An  aspirated,  teflon-coated  sensor  has  a  time  constant  of  about  30  minutes, 
but  this  varies  from  sensor  to  sensor.  Therefore,  the  humidity  differ¬ 
ence  read  between  two  sensors  will  be  unreliable  If  the  humidity  Is 
changing  faster  than  about  lOZ  RH  per  hour.  A  more  serious  deficiency 
appears  If  the  sensors  are  exposed  to  humidities  above  95Z  RH.  At  those 
high  humidities,  the  sensors  become  sluggish  and  exposure  to  99Z  EH  may 
cause  the  sensors  to  be  useless  for  as  long  as  several  days.  Further 
uncertainties  are  Introduced  by  the  use  of  long  cables  in  the  shipboard 
system  and  by  the  accumulation  of  sea  salt  on  or  near  the  sensors. 
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TABLE  A-2  Hygrodynanlcs  Digital  ZZ  Calibrations 


Wet-Dry 

Bulb 

%RH 

Ave 

Hygro 

%RH 

1 

2 

3 

4 

5 

6 

7 

8 

33. 

31.6 

1.4 

-1.4 

-0.6 

1.5 

-0.2 

-1.6 

-0.9 

1.8 

36. 

37.2 

-0.5 

-0.8 

-0.3 

1.6 

0.5 

-1.4 

-0.1 

1.2 

44. 

41.2 

-0.8 

-0.4 

-0.3 

0.6 

0.9 

0.4 

-0.1 

0.3 

NA 

41.5 

-0.5 

-0.4 

-0.3 

0.5 

0.9 

-0.4 

-0.1 

-0.1 

NA 

42.6 

-1.2 

-0.5 

-0.5 

-0.9 

1.8 

-0.4 

-0.2 

0.1 

NA 

65.8 

-0.2 

-1.4 

-2.5 

0.8 

1.8 

-0.3 

1.3 

0.3 

74.8 

73.6 

0.2 

0.2 

-1.3 

0.5 

-0.2 

0.4 

-0.4 

0.8 

88.3 

85.0 

-0.7 

0.9 

-0.4 

0.1 

-0.6 

0.0 

0.4 

0.2 

82.6 

86.7 

-1.4 

-0.6 

-0.5 

0.1 

-0.2 

0.3 

0.7 

0.3 

82.6 

87.0 

-1.0 

0.5 

-0.7 

0.0 

0.3 

0.7 

0.0 

0.2 

85.0 

88.0 

-0.4 

1.1 

-0.1 

0.1 

-1.8 

1.1 

-0.6 

-0.1 

93.9 

89.6 

0.1 

0.8 

-0.2 

0.3 

-2.2 

2.1 

-1.5 

0.4 

94.4 

92.3 

-0.6 

1.1 

-0.5 

-0.3 

-0.4 

1.2 

-0.5 

0.4 

94.4 

93.6 

-0.8 

1.0 

-0.8 

-0.3 

0.4 

-0.2 

0.8 

-0.4 

Ave  Deviat. 

-0.5 

-0.0 

-0.6 

0.3 

0.1 

0.1 

-0.1 

0.4 

0 

0.7 

0.9 

0.6 

0.7 

1.1 

1.0 

0.7 

0.6 

The  average  deviation  of  the  wet-dry  bulb  from  the  Ave  Hygro  reading 
0.4%  RH  with  a  standard  deviation  of  2.9%  RH. 
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3.  Velocity 

The  mean  wind  velocity  sensors  (Thomthwaite  cup  anemometers) 
were  calibrated  in  a  wind  tunnel  with  a  nominal  one  meter  square  test 
section.  The  cups  were  calibrated  four  at  a  time  with  a  symmetrical 
placement  to  reduce  relative  velocity  differences.  A  summary  of  the 
calibration  is  shown  in  Figure  A-6  expressed  as  individual  cup  velocity 
deviations  from  the  wind  tunnel  velocity.  From  this  data  it  can  be 
concluded  that,  for  gradient  purposes,  the  cup  anemometers  were  accurate 
to  about  i  42.  The  deviation  of  the  calibration  curve  at  low  wind 
speeds  is  due  to  errors  in  the  determination  of  wind  tunnel  velocity 
values . 

B.  Fluctuation  System 

The  essence  of  the  fluct\iatlon  calibrations  is  contained  in 
establishing  the  relation  between  the  instrumental  voltage  variation,  V, 
about  mean  voltage,  V.  with  respect  to  the  atmospheric  fluctuation,  X', 
about  its  mean  value  X,  i.e. 


X'  -  F  (  V,  X  )  V 


(1) 


The  cslibratlon  factor  F  (  V,  X  )  may  be  a  function  of  individual 
sensors  as  well  as  the  Instrument  being  used.  This  process  is  usually 
accomplished  by  measuring  the  dependence  of  X  on  V  and  relating  it  to 
F  (  V,  X  )  through  the  derivative. 


F  (  V,  X  ) 


(2) 
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I.  Temperature 

The  appropriate  temperature  fluctuation  expression  la: 


T'  -  (  a  G  R  V'  (3) 

where  a  Is  the  temperature-resistivity  coefficient  of  the  platinum  sensor 
wire,  G  Is  the  gain  factor  of  the  thermosonde  (set  at  G  ■  20  volt/(2),  and 
R  Is  the  resistance  of  the  sensor.  In  a  previous  work  (Schacher  and 
Falrall,  1976)  It  was  found  that  a  .0036  ^  with  little  variation 

from  sensor  to  sensor.  The  frequency  response  of  the  thermosonde  units 
was  about  1  kHz. 

2.  Humidity 

The  Lyman-alpha  sensor  measures  humidity  as  a  function  of  the 

0 

absorption  of  ultraviolet  light  at  the  1215A  Lyman-alpha  transition  of 
hydrogen  in  water  vapor.  It  consists  simply  of  a  UV  source  tube  and 
detector,  separated  by  an  absorbing  air  gap. 

The  detector  voltage  is  proportional  to  the  total  light  trans¬ 
mitted  and  Is  related  to  the  absolute  humidity,  q  (in  mb),  by  the 
equation 

V  •  exp  (-Xq)  (4) 


The  Lyman-alpha  Is  a  fast  response  device  and  Is  specifically  used  for 
measurement  of  humidity  fluctuations,  q' 


r 

V 


(5) 


The  quantity  in  Eq.  (4)  is  the  dry-air  voltage  which  can  be  found  by 
bathing  the  detector  in  dehumidified  air.  Since  water  soluble  windows 
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(LIF  or  MgF^)  are  required  for  transmission  in  the  ultraviolet,  is  a 
highly  variable  quantity,  a  property  which  makes  the  Lyman^alpha  a  poor 
device  for  measurement  of  total  humidity,  q.  However,  Eq.  (5)  shows  that 
the  fluctuating  part,  q',  depends  only  on  X,  which  is  a  property  of  water 
vapor  and  the  gap  setting  only  and  is  independent  of  the  window  condition 
so  that  q'  can  be  reliably  measured.  The  Lyman-alpha  was  calibrated  in 
dry  air,  ambient  air,  and  in  the  humidity  chamber  and  it  was  found  that 
X  ■  .225  mb~^  for  source  -  50  ya  and  a  1  cm  air  gap  (Fig.  A-7), 

Feasibility  of  Che  Lyman-alpha  sensor  for  shipboard  humidity  fluc¬ 
tuation  measurments  is  subject  to  considerations  of  the  survivability 
of  Che  windows  in  Che  ocean  environment  of  high  humidity,  sea  spray,  and 
rain.  In  Che  laboratory  humidity  chamber,  was  found  to  be  reduced  by 
a  factor  of  four  after  a  three-day  exposure  to  85Z  humidity.  Consequently, 
it  was  necessary  to  provide  physical  protection  for  the  sensors  during 
routine  operations. 

3.  Velocity 

The  hot  wires  (or  hot  films)  were  calibrated  using  a  TSl  1125 
calibrator.  It  was  Important  to  use  the  actual  cables  from  the  shipboard 
system  (or  their  equivalent)  to  prevent  errors  due  to  voltage  drops  in 
Che  cable.  The  appropriate  mean  dependence  function  is 


V  (U)  -  (V  ^  +  B  /T"  )**  (6) 

o 

2 

where  and  B  are  obtained  for  each  sensor  at  a  specific  overheat.  The 
derivative  of  (6)  yields  the  fluctuation  calibration  factor 


U' 


(7) 
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The  frequency  response  of  Che  velocity  bridge  was  about  20  kHz.  The 
sensor  frequency  response  depends  on  the  sensor  diameter  and  overheat, 
but  was  typically  greater  than  1  kHz  (Falrall  and  Schacher,  1977). 


V.  Data  Acquis It jog 

Three  criteria  applied  when  deciding  the  methods  for  acquiring  data: 

1)  data  was  to  be  taken  for  a  wide  range  of  local  meteorological  con¬ 
ditions,  2)  a  given  condition  should  be  tested  several  times,  preferably 
on  different  days  and/or  locations  to  test  the  statistical  validity  of  the 
data,  and,  most  Importantly,  3)  shipboard  conditions  must  be  such  that 
valid  data  can  be  obtained. 

The  most  Important  consideration  to  insure  validity  of  Che  data  was 
the  airflow  over  the  ship.  Since  all  of  the  sensors  were  located  close 
to  the  bow  of  the  ship  the  relative  wind  had  to  be  from  the  bow  or  the 
ship  could  severely  Influence  the  data  obtained.  It  was  found  that 
data  could  be  taken  when  the  wind  was  within  30*^  of  the  bow.  When  the  wind 
was  more  than  AS°  from  alignment  with  the  bow  the  wind  profiles  were 
noticeably  affected.  The  position  of  the  sun  was  a  second  consideration. 
When  the  sun  was  low  on  the  horizon  and  behind  the  ship  the  aspirators 
for  levels  1  and  2  were  in  shadow  and  those  for  levels  3  and  4  were 
In  sunlight.  Under  this  condition  the  temperatures  at  the  upper  levels 
were  elevated  enough  to  make  the  obtained  temperature  profiles  Invalid. 

Since  data  was  acquired  on  a  fluctuating  system,  a  time  average  had 
to  be  performed  to  obtain  final  parameters.  Averaging  over  10-minute 
and  20-mlnuce  time  spans  was  accomplished  depending  on  the  situation. 

Thus,  data  had  to  be  acquired  In  such  a  way  that  the  conditions  remained 
constant  over  a  time  at  least  as  long  as  the  planned  average  time.  This 
means  that  there  had  to  be  close  coordination  between  scientific  personnel 
and  ship  crew  so  that  no  course  changes  occurred  during  a  data  taking  run. 
Of  course,  changes  In  local  conditions  also  affected  the  data,  such  as 
moving  through  the  edge  of  a  fog  bank.  Even  seemingly  Inconsequential 
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occurrences  such  as  moving  Into  Che  lee  of  an  Island  when  the  ship  was 
several  miles  offshore  could  affect  the  data.  A  running  log  of  all 
ship  maneuvers  and  position,  and  of  local  conditions  was  maintained  and 
no  averages  were  performed  for  which  it  was  suspected  Chat  any  condition 
changed  during  Che  time  period  under  consideration. 

Several  tests  were  made  and  it  was  found  that  data  could  be  taken 
both  with  Che  ship  at  rest  and  underway.  However,  some  of  the  subsequent 
calculations  then  required  chat  the  true  wind  over  the  surface  of  the 
water  be  known  accurately.  Uncertainties  as  to  the  ship's  actual  speed 
make  calculation  of  true  wind  from  Che  relative  wind  and  ship  heading 
and  speed  too  inaccurate.  Thus  Che  ship  was  stopped  frequently  in  order 
to  obtain  Che  true  wind. 
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VI.  Data  Analysis 


The  analysis  procedures  changed  considerably  as  the  Instrumentation 
evolved  towards  Increased  automation.  Averaging  of  mean  data,  initially, 
a  painstaking  hand  process.  Is  now  performed  by  MIDAS  at  specified  time 
Intervals  (typically  10,  20  or  30  minutes) .  Rather  than  describing  all 
techniques  that  have  been  used,  discussion  Is  limited  to  those  of  more 
recent  application. 

A.  Mean  Profile  Analysis 

The  mean  profiles  and  gradients  were  established  under  the  assumption 
of  the  logarithmic  height  dependence. 


X  *  A  log  Z  +  B 


(8) 


The  gradient  at  the  height  Z  is 


dz 


A 

Z 


(9) 


In  the  case  of  temperature,  one  Is  Interested  In  the  potential  temperature, 

e. 

F  -  T  +  .0098  Z  (10) 


and  the  virtual  potential  temperature,  9^, 

9  -  e"  +  .  61  q  T 

V  ^ 


(11) 


where  q  is  the  specific  humidity  in  grams  of  water  vapor  per  gram  of  dry 
air.  The  actual  fits  of  the  data  were  done  either  by  the  standard  least 
squares  method  or  by  subjective  ("eyeball")  graphical  techniques.  The 
presence  of  Individual  erroneous  values  due  to  printer  errors,  mlcropro- 
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cessor  errors,  and  Instrumental  glitches  required  the  data  to  be  subjec¬ 
tively  edited  even  in  the  least  squares  analysis.  In  either  case,  one  is 
required  to  identify  and  reject  "bad"  values  -  a  delicate  and  unsatisfying 
process.  Figure  A-8  shows  examples  of  profiles  in  temperature,  humidity  and 
velocity. 


B.  Turbulence  Analysis 

There  are  four  different  methods  available  to  calculate  turbulence 
parameters  from  the  fluctuations  in  the  atmospheric  variables:  spectral, 
derivative  (dissipation),  difference,  and  RMS.  The  spectral  method  is 
based  upon  the  assumption  of  "local  isotropy"  and  the  Kolmokorov  -3/3 
slope  of  the  one-dimensional  power  spectral  density,  , 

<|)^(k)  -  .25  k  (12) 


where  C  is  the  structure  function  for  the  parameter  x  and  k  is  the  wave 
number.  Using  Taylor's  "frozen  turbulence"  hypothesis  (k  *  2iTf/U), 
is  found  by  performing  a  fourier  spectrum  analysis  of  a  signal  in  the  fre¬ 
quency  domain  (f) , 

-2/3 

f  <t>^  (f)  “  k  (k)  -  .25  C^“  I  (13) 


Therefore, 

c/  -  (f)]  (14) 

which  is  expected  to  apply  in  the  frequency  range  from  about  0.1  to  100  Hz. 
(This  is  the  inertial  subrange  of  turbulence.)  The  structure  function  can 
also  be  found  by  measuring  the  variance  of  the  difference  in  the  variable 
X  at  two  points  separated  by  a  known  distance,  d, 

-  [X  (r)  -  X  (r  +  d)]^  (15) 
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Defined  in  this  manner,  is  Independent  of  d  in  the  Inertial  subrange. 

The  RMS  method  is  based  upon  measuring  the  variance  of  the  signal  fluc- 

2 

tuations  between  selected  frequency  limits  and  relating  this  to  through 
Eq.  (13): 


k 

/  “ 


(k)  d  k 

X 


r.2 


(X’  )■ 

rms 


(16) 


Again  assuming  k  *  2iTf/U,  then 


8  /  2ir 


3  I  U 


2/3  (X’ 


rms 


)^ 


(f  -  f^-/^) 


(17) 


where  f  and  f  are  the  upper  and  lower  frequency  limits,  and  X' 

u  c  rms 

represents  the  measured  RMS  fluctuations  of  the  filtered  signal.  Since 

Che  rate  of  velocity  turbulence  dissipation,  £,  is  traditionally  used, 

2 

rather  chan  C  ,  one  should  use  the  relation 
u 

2  2/3 

*  2.0  (18) 


to  calculate  e.  Due  to  the  extreme  sensitivity  to  high  frequency  noise,  it 

was  found  chat  the  derivative  method  was  unsatisfactory  for  shipboard  data. 

2  2 

The  spectral  method  was  used  to  calculate  C  and  e  (data  shewn 

T  q 

in  Figure  A-9).  This  method  has  the  advantage  that  the  quality  of  the  data 
is  immediately  obvious  and  noise  spikes  (60,  120,  130  Hz)  can  be  ignored. 

The  primary  disadvantage  of  the  spectral  method  is  the  time  required  to 
analyze  four  levels  of  data  with  two  or  three  signals  for  each  level.  Al¬ 
though  it  requires  two  microthermal  sensors  per  level,  the  difference  method 

2 

has  proven  to  be  an  Ideal  teennique  for  measuring  because  the  RMS  temp¬ 


erature  difference  fluctuations  can  be  automatically  averaged  using 


Che  ADC  capability  of  MIDAS.  Since  this  is  a  temperature  difference 
measurement,  it  is  less  subject  to  external  sources  or  error  (such  as 
humltlty  -  salt  induced  temperature  spikes  (Frlehe,  1977)).  The  RMS 
technique  has  proven  to  be  quite  satisfactory  for  measuring  e,  since  it 
too  lends  Itself  to  automatic  averaging  on  MIDAS.  Comparison  of  spectral 
and  RMS  measurements  of  e,  through  the  quantity  ■  (K  Z 
(Fig.  A-10)  are  quite  consistent.  is  the  friction  velocity,  K  is 
Von  Karmon’s  constant  (.35)  and  Z  is  Che  height  above  the  surface.  This 
expression  is  valid  for  near  neutral  stability  where  most  of  the  data 
was  taken. 
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VII.  Experiments  Performed 

Observational  experiments  were  conducted  aboard  the  R/V  Acanla 
(Figure  A’-l)  during  cruises  to  San  Nicholas  Island  (SNI),  in  the  vicinity 
of  SNI,  In  Monterey  Bay,  and  on  an  open  ocean  cruise  in  the  fall  of  1976. 
Dates  and  locations  of  the  experiments  which  yielded  results  presented 
here  are  listed  in  Table  A-3. 

The  route  taken  to  and  from  SNI  appears  in  Figure  A-11  for  that  portion 
of  the  cruise  for  which  data  was  taken.  The  ship's  anchorage,  position, 
and  track  In  Che  vicinity  of  the  northwest  tip  of  SNI  appear  In  Figure  A-12. 
For  Monterey  Bay  results,  all  data  were  collected  while  Che  R/V  Acanla  was 
anchored.  The  general  anchorage  locations  which  were  used  in  Monterey 
Bay  appear  In  Figure  A-13.  The  anchorage  locations  were  changed  to  accom¬ 
modate  coincident  ship  to  shore  or  shore  to  shore  optical  propagation 
experiments.  Wind  directions  during  the  Monterey  Bay  experiments  were 
generally  from  the  north  to  northwest  so  the  location  of  anchorage  is  not 
expected  to  be  a  factor  in  the  results  and  they  are  not  distinguished  with 
regard  to  location.  Details  pertaining  to  the  1976  open  ocean  cruise  appear 
in  a  separate  report,  Fairall,  et  al.  (1978). 

Changes  were  made  in  Che  shipboard  measurement  systems  and  sensor 
mounting  arrangements  during  the  sequence  of  experiments  associated  with 
these  results.  In  general,  mounting  arrangements  were  changed  to  improve 
the  vertical  resolution  of  the  measured  parameters.  During  the  first  SNI 
cruise  (1-4  March  1973),  fluctuation  (temperature  and  velocity)  measure¬ 
ments  were  made  at  only  two  levels  and  mean  wind,  temperature  and  humidity 
measurements  were  made  only  at  one  level.  Sea  surface  temperature  was  also 
measured.  During  the  second  SNI  cruise  (19-23  September  1973),  both  fluc¬ 
tuation  and  mean  measurements  were  made  at  four  levels  which  differed 
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slightly  from  the  final  configuration  shown  In  Figure  A-1.  All  the  Honterey 
Bay  experiments  were  performed  with  the  final  mounting  configurations 
described  above.  Sensors  described  for  fluctuation  and  mean  measurements 
were  used  In  all  experiments  but  the  recording  systems,  particularly  that 
for  the  mean  data,  were  changed. 

Additional  discussions  on  analyses,  instrumentation  and  experimental 
conditions,  along  with  additional  results,  are  contained  In  theses  listed 
In  the  references,  which  utilized  data  from  these  experiments. 


! 
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TABLE  A-3 


SUMMARY  OF  OBSERVATIONAL  EXPERIMENTS  TO 


Date 

Location 

1973 

1-5  Mar 

SNI 

18-21  Sep 

SNl 

15-19  Oct 

Mtry  Bay 

26  Nov-5  Dec 

Mtry  Bay 

1974 

15-18  Jan 

Mtry  Bay 

25-26  Feb 

Open  Ocean 

27-29  Feb 

Mtry  Bay 

6-7  Mar 

Mtry  Bay 

13-15  Mar 

Mtry  Bay 

25-28  Mar 

Mtry  Bay 

17-21  Jun 

SNI 

12-16  Aug 

Mtry  Bay 

16-21  Sep 

Open  Ocean 

18-22  Nov 

Mtry  Bay 

1975 

9-10  Jan 

Mtry  Bay 

24-28  Mar 

Mtry  Bay 

20-24  May 

Mtry  Bay 

25-29  Aug 

Mtry  Bay 

11-16  Sep 

Open  Ocean 

9-11  Dec 

Mtry  Bay 

1976 

29  Mar-2  Apr 

Mtry  Bay 

26-30  Apr 

Mtry  Bay 

26-30  Jul 

Mtry  Bay 

20  Sep-14  Oct 

Open  Ocean 

1977 

15-18  Feb 

Mtry  Bay 

FEBRUARY  1977 

Prevailing  Condition 

Unstable 

Neutral-Stable 

Unstable 

Unstable 

Unstable 

Stable 

Stable 

Stable 

Stable 

Unstable 

Neutral 

Neutral -Uns  tab le 
Unstable 

Neu  tr al-Uns  tab 1 e 

Stable 

Unstable 

Neutral 

Neutral 

Unstable 

Unstable 

Unstable 

Neutral -Uns tab le 
Neutral 

Neutral-Unstab le 
Neutral 
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Vlll.  System  Evaluation 


A.  Mean  System 

Accurate  measurements  of  atmospheric  surface  gradients,  especially 
In  the  marine  boundary  layer,  are  extremely  difficult.  The  severity  of 
the  problem  can  be  Illustrated  by  noting  typical  conditions  In  the  marine 
boundary  layer  (Table  A-A)  in  terms  of  the  typical  gradients  OX/3Z 
Z  *  10  meters),  scaling  parameters  (X^)  and  total  difference  from  Z  -  5  to 
Z  *  15  meters  (AX  >  |X  (IS)  -  X  (5)  |).  The  scaling  parameter  Is  defined  by 


ax/3z  .  ^  f^ 


(19) 


where  f ^  is  a  stability  corrections  factor  (gusinger,  et  al.,  1971).  To 
measure  X^  to  lOZ  accuracy  requires  the  single  level  accuracies  Indicated 
In  the  table,  compared  to  our  present  accuracies.  The  temperature,  T,  Is 
used  In  the  table  rather  than  the  virtual  temperature  because  It  is  the 
variable  actually  measured. 

Although  the  gradients  are  not  required  to  calculate  the  turbulent 

parameters,  they  are  of  utmost  importance  In  evaluating  the  theoretical 

behavior  of  the  boundary  layer  through  specification  of  the  stability 

(Richardson  number)  and  through  normalization  of  the  turbulence  parameter, 
2 

(Wyngaard,  et  al.,  1971).  This  lack  of  gradient  accuracy  led  to  con¬ 
struction  of  a  near  surface  buoy  system  to  exploit  the  Z~^  dependence  of 
the  gradients.  With  the  present  systems,  the  gradient  measurements  repre¬ 
sent  the  weak  link  In  efforts  to  evaluate  the  turbulent  boundary  layer 
similarity  expressions  In  the  marine  environment. 


B.  Fluctuation  System 

Despite  a  large  array  of  factors  Influencing  the  performance  of  the 
turbulence  measurement  system  (noise,  ship  effects,  salt  on  sensors. 


TABLE  A-4 


Typical  Marine  Boundary  Layer  Profile  Parameters 


T,  ®C 

0,  Vsec 

q,  ginAgm 

met 

-  .020 

.040 

-  .025 

-  .07 

.15 

-  .08 

AX 

.20 

10  % 

5  % 

Present  Accuracy 

? 

5  % 

3  % 

Required  Accuracy 

.04 

2  % 

1  % 
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calibration  accuracy)  the  values  of  turbulence  parameters  were  quite  con- 
slstent.  Figure  A-14  is  a  time  history  of  four  levels  of  and  E  showing 

the  expected  height  dependence.  Figure  A-15  is  a  similar  plot  for  three 
levels  of  e.  The  scatter  Inherent  in  the  e  data  is  nicely  illustrated 
by  removing  the  height  dependence  in  a  time  history  (Fig.  A-15)  of  buoy 
data.  Note  that  the  single  level  scatter  is  much  smaller  than  the  level- 
to-level  disagreements.  The  average  profile  estimate  (based  on  Eq.  19) 
of  is  about  .12  m/sec  whereas  the  average  turbulence  estimate  is  about 
.10  m/sec.  This  is,  in  fact,  good  agreement  and  Illustrates  the  superior¬ 
ity  of  the  buoy  systea  for  gradient  measurements.  The  compatabillty  of 
gradient  and  turbulence  measurements  can  be  examined  in  calculations  of 
normalized  heat  flux,  F^,  using  the  profile  method 

(30/32)  (20) 

and  the  turbulence  method 

Fjj  -  ±  e  F  (Rj^)  (21) 

where  K,^  and  F(R^)  are  functions  of  stability.  These  methods  have  been 
favorably  compared  (Falrall,  et  al.,  1978)  for  fog  events  during  a  marine 
fog  study  called  CEWCOM-76  (Fig.  A-16). 


C.  Summary 

It  is  apparent  that  this  line  of  research  is  permeated  with 
Instrumental  difficulties.  The  atmospheric  gradients  over  the  ocean  are 
small  enough  so  that  even  state-of-the-art  mean  measurements  give  fairly 
poor  accuracy.  The  turbulence  measurements  are  Inherently  well  above 
noise  levels  but  the  sensors  are  subject  to  an  array  of  environmental 
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problems:  salt  and  water  films,  platform  motion,  aerodynamic  flow  dis¬ 
tortion  by  the  ship,  radio  and  radar  interference,  and  even  corrosion. 
Despite  these  problems,  the  syst^  described  In  this  paper  has  yielded 
valuable  and  unique  data  about  the  open  ocean  boundary  layer.  Those 
theoretical  discrepancies  that  still  exist  will  hopefully  be  resolved 
through  the  continual  evolution  and  improvement  of  Instrumentation  and 
measurement  techniques. 
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APPENDIX  B 

MOS  STABILITY  FUNCTIONS 


The  forms  of  the  mean  gradient  functions  (Businger  et  al.,  1971). 


♦u(0  -  (1  -  155)'''''* 

c<o 

*  (1  +  4.75) 

C>0 

♦^(5)  «  (1  -  195)“''^^ 

5<0 

♦^(5)  *  (1  +  6.45) 

C>0 

The  mean  profile  function 

'i>U)  =  l-<^(5)  -  31n  ^(5)  +  21n  +  2 

2 

tan’^  ^U)  -  ir/2  +  In  (Itiliil; 

The  dimensionless  velocity  dissipation  function 

(Wyngaard  and  Cote,  1971) 

Ey(5)  *  (1  +  .5l5^/^)^/2 

C<0 

E^J(5)  »  (1  +  2.55  2/3)3/2 

OO 

The  dimensionless  temperature  structure  function  parameter  (Wyngaard  et  al, 
1971). 

f^(5)  =  4.9(1  -  75  )’^^^ 

5<0 

fjU)  *  4.9(1  +  2.45 


C>0 
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TABLE  C-1,  EM  Propagation  data  (bulk  calculations). 
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TABLE  C-2  MOS  scaling  parameters  and  EM  propagation  results.  Comparl 
son  of  the  bulk  method  (first  line)  with  turbulence  U«  and 
profile  T*,  q*  (second  line). 


Z/L  Monin-Obukhov  stability  parameter,  c 
Velocity  scaling  parameter  (m/sec) 

Potential  temperature  scaling  parameter  (°C) 

Water  vapor  mixing  ratio  scaling  parameter  (gm/kg) 
Duct  height  (m) 

Fm  Minimum  frequency  trapped  (GHz) 
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18.3 

13.0 


5 


3/0  5 


5 


■1.98E-02 
>1.  38E-01 


0.35 

0.46 


•0.03 

•0.18 


■0.18 

■0.33 


17.7 

20.0 


-5.02E-02 

0.33 

-0.05 

-0.21 

17.5 

4.-0 

170 

3/0  5 

2400 

-1. 32E-01 

0.46 

-0.17 

-0.33 

20.0 

4.0 

171 

3/06 

30 

-8. 255-02 

-7. 35E-02 

0.33 

0.50 

-0.10 

-0.15 

-0.27 

-0.  30 

19.3 

21.0 

•  •• 

3.0 

-2.93E-02 

0.34 

-0.05 

-0.  19 

17.9 

4.7 

172 

3/06 

100 

-7.03E-02 

0.60 

-0.17 

-0.31 

21.3 

-2.495-02 

0.34 

-0.04 

-0.21 

1^.7 

4.1 

173 

3/0  5 

130 

-1. 345-01 

0.4  7 

-0.17 

-0.  32 

10. 

4 . 1 

-3.585-02 

0.34 

-0.0  3 

-0.22 

19.0 

4.1 

174 

3/06 

200 

-6.22E-02 

0.58 

-0.17 

-C.31 

22.4 

3.4 

6 


-3.28E-02 


0.32 


.2 


0.03 


-0.28 


24.6 


#  Date  GfiT 


21% _ 3/LQ _ L810_ 


L76. _ 3/ijQL _ laQQ-- 


_2i7 _ 3/10 _ laaa. 


.278 _ a/lQ. _ 10  OQ. 


_ 27J> _ 3/X.O.  . .  J103(L _ 


L5^6ilEr.Q2_ 

■6.13E-03 


0.31 


r(l..Q3- 

0.00 


-2.40E-02  0.31  -0.00 


l.4^33E=Q2_ 

5.31E-01 

a^4.4£.--Q2- 

•1.45E-02 

r2..19E=L02_ 

•2.28E-02 


.JI..24 _ 

0.32 

_0.27. _ 

0.40 

0.32 


-3.75E-02  0.31 


^.01- 

0.24 

l(U.OQ_ 

0.00 

r.0.00 

0.00 

•n.on 

-0.00 


-0.04 

_.^o..a4__ 

-0.07 

0.17 

-0.07 

...--0^2_ 

-0.07 

__=Q.J13_ 

-0.07 


281  3/10 


-1.16E-02  0.36  0.01  -0.06 


.3LJ. 

29.4 


7.3  18.4 


-43^0 _ 

14.3 

_ 

15.4 


7.4  18.0 


6.7  20.3 


-5,35E-02  0.29  C.Ol  -0.15  13.5  7.: 

-6.525-02  Q.  23  -0.00  -0.15  .  1.1^4  7. 

-1.12E-G1  0.29  -n.Ol  -0.19  14.6  6.5 


ttIWi 


-3.26E-02  0.31  0.02  -0.20  19.2 


-4.68E-02  0.31  0.02  -0.20  17.9 


-7.26E-02  0.30  0.00  -0.20  16.2  5.5 


9.72E-03  0.30  0.03  -0.13  18.2 


6 v42E-03  -  Ok31  0,03  -0. 17  23.2 


WkVHHBc  r4  ■  r  ff  1 1 11 


3.45E-02  0.29  0.04  -0.12  22.6 


6 . 69B-02  0 , 29  0 .06  -0. 15  137, 9 


7.07r-02  0.30  0.07  -0.15  140,5 


1.27E-01  0,30  0.09  -0,10  709,1 


-1.033  01 


ritrn^Tiwi 


-6.53  -1.25 


1,17E-01  0,35  0.12  -0,14  543,3 


Ei—H 


6,40E-02  0.38  0,07  -0.14  74.1 


4.37E-02  0.30  0.05  -0,16  48,9 


2,7lE-02  0.30 


i 


5  -0.17  34.8 


-1.04E-02  0,28  0,03  -0.22  24.1 


E»1 


■1.55S-01  0.30  -0,03  -0.34  21.3 


■1.16B-01  ^  0,32  -0,03  -0.38  23.3  3.2 

.2fi  15.6 5.9 


.32  20.9  3.8 


•1,52E-01  0.30  -0.03  -0.32  20.0  - 


■1.90E-01  0.32  -0.05  -0.43  24.2 


•1.13E-01  0.41  -0.07  -0.43  26.7  2.$ 


VlHiE  I M I  r  J 


fcTiTil 


•2.44E-01  0.37  -0.08  -0.69  32.6 


•3.67B-01  0.36  -0.12  -1.09  42.0  l.J 


4.87E-01  0.30  -0.12  -1,31  45.2 


4.99BH)1 


:  1  jCri  1 I  iv  1 1 


IH 


-0.12  -1,56  51.1  1,01 


4.77E-01  0.30  -0.15  -1.31  45. 


4.24E-01  Q.31  -O.Il  -1.28  45.9 


6.29E-01  0.20  -0.19  -1.36  43.3 


iilvrJ 


-4.90E-0t  0.20  -0.12  -1,26  41.6 


-5.05E-01 


0.22  -0.01  -0,28 


14.0 


6.9 


TABLE  C- 


3,  Similar  to  TABLE  C-2,  but  for  Bouy  Data  Only. 
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1.74E:-01 


9.  20E-02 

i.aiE-oi 


3.32E-02 

1.84E-01 


1.06E-01 

1.30E-01 


3/11 

3/11 


10  30 

1100 

1130 


■8.86E-02 

■4.68E-02 

•8, 12E-02 
•7.25E-02 

■8.04E-Q2 

■7.72E-02 


0.23 

0.31 

0.24 

0.30 

0.24 

0.30 


•0.01 

0.02 

•0.01 

0.00 

•0.01 

0.00 


lO-jJLL 

-0.20 

•0. 16 

■n.  20 

•0.17 

•0.22 


14.2 
17.9 

13.4 

16.2 

13.7 

17.6 


292 

3/11 

12  30 

-7.47E-02 

0.  26 

-0.01 

-0. 16 

13.5 

7.3 

■  ; 

-2.31E-02 

0.31 

0.02 

-0.  18 

18.  5 

4.5 

4 

293 

3/11 

1300 

-1. 02E-01 

0.25 

-0.01 

-0.  25 

18.3 

4.6 

V" 

;4 

9.72E-03 

0.30 

0,03 

-0.13 

18.2 

4.6 

1 

u  ^ 

1 

Date 

CMT 

Z/L 

u* 

T* 

a* 

Z* 

Ft 

i  294 

3/11 

1330 

-8.65E-02 

0.25 

-0.00 

-0.25 

19.2 

4.3 

6.42E-03 

0.31 

0.03 

-0.  17 

23.2 

3.2 

295 

3/11 

1400 

-1.62E-01 

0.22 

-0.0  2 

-0.27 

-  17^.. 

4.9 

3.45E-02 

0.29 

0.04 

-0.12 

22.6 

3.4 

1  296 

3/11 

1430 

-1.56E-01 

0.23 

-0.0  3 

-0.  24 

15.7 

5.8 

1 

1 

6. 53E-02 

0.29 

0.06 

-0.15 

137.0 

0.2 

'1 

'  2  97 

3/11 

1500 

-1.33E-01. 

-O.Od 

—0 .23 

_ 

5^6 _ 

7.07E-02 

0.30 

0.07 

-0.  15 

148.5 

0.2 

208 

3/11 

1530 

-1.24E-01 

0.26 

-0.02 

-0.28  . 

19.  3. 

4.J3 _ 

1.  27E-01 

0.30 

0.09 

-0.  10 

70^.  1 

0.0 

299 

_3/ll 

1600 

-.uoi^ai _ 

0.28 

-0.02 

-0.27 

_  19.  1 

4.3 

!  '  i 

-1.03E  01 

0.31 

-6.53 

-1.25 

2.2 

111.2 

300 

3/11 

1700 

-2.43E-02 

0.31 

..  _D.02 

-0.22 

22.1 

_ _ 

•  t 

1.17E-01 

0.35 

0.12 

-0. 14 

543.3 

0.0 

301 

3/11 _ 

1730 

_7  92p._02 

0-30 

-Q.Ql  . 

-0-27 

20.4 

3-9 

6.40E-02 

0.38 

0.07 

-0.14 

74.1 

0.6 

302 

3/11 

.mo _ 

CL..26 

-0.01 

_-0..27 _ 

JLSUl 

4.3 

4.37B-02 

0.30 

0.05 

-0. 16 

48.9 

1.1 

i  303 

J/11 _ 

1830 

6  S'"  01  - 

.(U27 

rO_iL.QA_._ 

-0...1S _ 

18.3 

4.6  . 

i 

2.71E-02 

0.30 

0.05 

-0.17 

34.3 

1.3 

i  304 

3/11 

1900 

-1. 942-01 

0.24 

-0.05 

-0.29.  . 

17.4 

5.0 

i 

! 

-1.04E-02 

0. 23 

0.03 

-0.  22 

24.1 

3.0 

;  305 

3/11 

1930 

-2.57E-0L 

0.23 

-0.07 

-C.29 

16-0 

.  5.6 

‘! 

-2.30S-01 

0.33 

-0.03 

-0.4  2 

22.9 

3.  3 

Date 


3/11  2200 


-2.52E-01 


-1.52E-01 


0.27 


0.30 


-0.11 


0.03 


-0.25 


-0.32 


14.3 


20.0 


3/11  2230 


-2.38E-01 

-l.'9dE-b'l 


0.28 

0.32' 


-0.05 


-0.43' 


15.9 

24.2 


3/11  2300 


-1.14E-01 

-1.13e-0l' 


0.40 

0.41' 


-0.10 

-0 . 0'7 


-0.43 


17.8 

26.7 


3/15  1130 


-2.01E  00 
-T.D8E  00 


0.11 

0.20 


-0.15 

-0.20 


zO-iA. 

■-0*.4  7 


12.7 

12.7 


3/15 


1230 


-4.62E-01 

-r7T5E-01 


0.21 

0.25 


-0.12 

-0.04 


-0.34 

-0:22 


15.6 
14.  1 


3/15  1300 


3/15 


3/15 


1400 


■1.590  00 


■5 . 0  5E-01 
■1.376  00 


. oOE— 0 
■1.70E  00 


•4.15E  00 


0. 12 


0.22 

0.11 


.21 

0.11 


0.21 


-0.14 


0.01 

-0.15 


0.03 

-0.14 


0.30 


-0.39 


-0.23 
-0.4  2 


-0.05 

-0.40 


-0.30 


13.0 


14.0 

13.1 


13.' 


17.5 


TABLE  C' 


•4,  Meteorological  data.  The  data  is  given  by  level  number 
(ie.,  U2  is  the  wind  speed  at  level  2). 

U  Relative  wind  speed  (m/sec) 

PHI  Relative  wind  direction  (degrees) 

TS  Sea  surface  temperature  (°C) 

T  Air  temperature  {°C) 

H  Relative  humidity  {%) 
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U3  U4  PHI  TS  T1  T2 

(n/sec)  (deg)  (centigrade) 
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